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Preface 


The  following  report  covers  a  discussion  of  the  research  done 
on  the  Mesosphere  Dynamics  Project,  Contract  No.  AF  1 9(604)-5492, 
during  the  life  of  the  project,  1  May  1959  to  30  November  1961. 

The  efforts  of  the  project  were  directed  towards  study  of  the 
dynamics  and  thermodynamics  of  the  stratosphere  and  mesosphere. 
Because  of  the  fundamental  role  played  by  ozone  to  problems  in  this 
region  of  the  atmosphere,  research  related  to  the  formation,  obser¬ 
vation,  and  distribution  of  ozone  became  of  central  interest  to  the 
project.  The  research  is  reported  in  three  parts; 

Part  I;  Accuracy  of  the  determination  of  the  vertical  distri¬ 
bution  of  atmospheric  ozone  by  the  infrared  method,  by  Katsuyuki 
Ooyama. 

The  vertical  distribution  of  atmospheric  ozone  as  determined 
by  the  infrared  method  is  limited  in  vertical  resolution  and  is  very 
sensitive  to  observational  errors.  A  set  of  linearized  error  equations 
was  derived  from  the  three  basic  equations  of  the  infrared  method, 
and  the  sensitivity  of  the  determined  ozone  distribution  to  errors  in 
the  observational  input  data  (total  ozone,  infrared  absorption,  and 
infrared  emission)  was  computed  for  a  few  cases.  It  was  shown  that 
in  the  region  above  the  tropopause,  the  ozone  distribution  is  most 
sensitive  to  measurements  of  infrared  absorption.  An  estimate 
was  made  of  the  accuracy  of  the  inverted  use  of  the  Curtis -Godson 
approximation  in  obtaining  the  mean  pressure  of  the  ozone  distri¬ 
bution. 


vii 


Part  II:  Absorption  processes  in  the  stratosphere  and  meso¬ 
sphere,  by  Julius  London  and  Cuddapah  Prabhakara. 

The  photochemical  equilibrium  distribution  of  ozone  was  com¬ 
puted  for  various  latitudes  for  summer  and  winter  and  the  absorption 
of  solar  radiation  by  ozone  and  molecular  oxygen  was  calculated.  It 
was  found  that  the  absorption  by  ozone  in  the  atmosphere  below  80  km 
was  everywhere  larger  than  that  by  molecular  oxygen.  The  region 
of  maximum  absorption  per  unit  mass  was  shown  to  coincide  with  the 
stratopause. 

Part  III;  The  distribution  of  total  ozone  over  the  Northern 
Hemisphere,  by  Julius  London. 

Maps  of  the  distribution  of  total  ozone  over  the  Northern 
Hemisphere  were  analyzed  for  the  whole  year  and  for  each  of  the 
four  seasons  from  data  accumulated  over  the  period  1926-1959.  The 
analysis  showed  marked  longitudinal  variations  of  total  ozone  even 
for  the  annual  average.  Ozone  ridges  appear  well  developed  over 
Eastern  North  America,  Eastern  Asia,  and  Central  Europe.  The 
hemispheric  distribution  of  total  ozone  was  shown  to  be  related  to 
the  atmospheric  circulation  in  the  upper  troposphere  and  lower 
stratosphere. 

In  addition  to  those  mentioned  above,  Mr.  James  Neiger 
and  Capt.  Robert  Simonson  contributed  subsi-antially  in  time  and 
effort  to  the  accomplishments  of  the  project. 
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PART  I 


Accuracy  of  the  Determination  of  the  Vertical  Distribution 
of  Atmospheric  Ozone  by  the  Infrared  Method 

by 

Katsuyuki  Ooyama 

ABSTRACT 

The  vertical  distribution  of  atmospheric  ozone  as  deter¬ 
mined  by  the  infrared  method  is  limited  in  vertical  resolution 
and  is  very  sensitive  to  observational  errors.  A  set  of  linear¬ 
ized  error  equations  was  derived  from  the  three  basic  equations 
of  the  infrared  method,  and  the  sensitivity  of  the  determined 
ozone  distribution  to  errors  in  the  observational  input  data 
(total  ozone,  infrared  absorption  and  infrared  emission)  was 
computed  for  a  few  cases.  It  was  shown  that  in  the  region 
above  the  tropopause,  the  ozone  distribution  is  most  sensitive 
to  measurements  of  infrared  absorption.  The  numerical  re¬ 
sults  indicate  that  a  positive  error  of,  e,  g.  ,  one  percent  in 
the  observation  of  infrared  absorption  could  erroneously  in¬ 
crease  the  ozone  amount  in  the  layer  tropopause  to  50  mb  by 
15  percent  and  erroneously  decrease  the  amount  of  ozone  in 
the  layer  50-30  mb  by  about  10  percent. 

In  the  infrared  method  the  Curtis -Godson  approximation 
is  used  in  an  inverted  way  in  order  to  obtain  the  mean  pressure  of 
the  ozone  distribution.  An  estimate  of  the  accuracy  of  the  in¬ 
verted  use  of  the  approximation  was  made  by  making  use  of  a 
simple  spectral  model,  and  the  effect  of  this  theoretical  error 
on  the  determined  ozone  distribution  was  also  included  in  the 
above  error  analysis.  The  mean  pressure  was  found  to  be  under¬ 
estimated  by  about  20  percent  in  the  case  of  atmospheric  ozone, 
and  the  resulting  ozone  distribution  would  be  severely  dis¬ 
torted  by  this  error  alone. 


1.  Introduction 


Although  the  total  amount  of  atmospheric  ozone  in  a  vertical 
column  is  now  observed  daily  at  many  stations  over  the  world,  the 
usefulness  of  ozone  observations  in  the  study  of  atmospheric  motion, 
especially  in  the  stratosphere,  would  be  radically  increased  if  the  ver¬ 
tical  distribution  of  ozone  were  observed  on  a  synoptic  basis.  Among 
the  several  methods  of  observation  which  have  been  developed,  direct 
methods  by  sending  apparatus  aloft  are  preferable  in  terms  of  poten¬ 
tially  obtainable  accuracy.  Nevertheless,  until  the  time  when  more 
and  inexpensive  ozonesondes  become  available,  indirect  methods  such 
as  the  Umkehr  and  the  infrared  methods  may  continue  to  be  important 
means  of  supplying  needed  data  of  the  vertical  ozone  distribution. 

The  infrared  method  as  discussed  in  this  paper  is  basically 
the  same  as  that  recommended  by  the  International  Ozone  Commission 
(1956)  and  as  described  by  Goody  and  Roach  (1956).  A  variation  of 
this  method,  particularly  of  the  computational  procedure,  was  pro¬ 
posed  by  Vigroux  (1959).  Empirical  comparison  of  the  vertical  ozone 
distributions  determined  by  Vigroux  with  those  which  were  simultaneous 
ly  determined  by  other  methods  --chemical  and  optical  ozonesondes, 
and  Umkehr  --were  recently  published  by  Brewer,  DQtsch,  Milford, 
Migeotte,  Paetzold,  Piscalar,  and  Vigroux  (I960).  A  visual  inspec¬ 
tion  of  the  results  reveals  that  the  values  of  the  ozone  density  at  any 
level  which  are  determined  by  the  different  methods  often  differ  by 
50  percent  or  more.  Goody  and  Roach  (1956)  showed  that  the  ozone  dis¬ 
tribution  in  the  stratosphere  as  determined  by  the  infrared  method  was 
extremely  sensitive  to  the  error  in  the  absorption  measurement  in  the 


9.6|j.  ozone  band.  As  a  result,  the  British  group  aimed  at  determination 
of  only  the  tropospheric  ozone  amount  from  the  measurement  of  emission 
intensity  in  the  9.6p  band  (Goody  and  Roach,  1958;  Walshaw,  I960). 

Another  method  of  determining  the  vertical  ozone  distribution 
from  infrared  observations  was  independently  developed,  and  a  long 
series  of  observations  at  Flagstaff,  Arizona  was  analyzed  by  Epstein, 
Osterberg,  and  Adel  (1956).  This  method  was  based  on  the  same  prin¬ 
ciples  as  that  of  Goody  and  Roach,  but  it  is  now  considered  obsolete 
(Epstein,  1957)  partly  because  of  an  improper  assumption  on  the  pressure 
dependency  of  the  infrared  absorption  by  atmospheric  ozone.  The  study 
reported  here  originated  in  our  effort  to  reanalyze  the  unique  set  of 
Flagstaff  observational  data  by  the  method  recommended  by  the  Inter¬ 
national  Ozone  Commission.  Unfortunately,  the  attempt  at  reanalysis 
had  to  be  abandoned  because  of  the  strong  suspicion  that  systematic 
errors  were  present  in  the  observational  data  for  which  no  reliable  cor¬ 
rection  procedure  could  be  found. 

The  purpose  of  this  paper  is  to  make  a  theoretical  analysis  of 
the  accuracy  of  the  determination  of  the  vertical  distribution  of  atmo¬ 
spheric  ozone  by  the  infrared  method.  The  error  analysis  of  the  deter¬ 
mined  ozone  distribution  is  made  in  terms  of  the  errors  in  observational 
data  and  also  in  the  inverted  Curtis -Godson  approximation.  Although 
technical  aspects  of  the  infrared  observations  are  not  discussed  here, 
the  result  of  the  error  analysis  will  be  a  guide  to  those  who  may  plan 
to  make  such  observations.  Since  the  ozone  distribution,  particularly 
in  the  stratosphere,  is  found  to  be  very  sensitive  to  observa.tional  errors, 
this  fact  should  be  borne  in  mind  by  those  who  are  mainly  interested  in 
results  obtained  by  the  infrared  method. 


2.  Formulation  of  the  infrared  method 


Absorption  due  to  the  9.6p.  band  of  ozone,  like  infrared  ab¬ 
sorption  due  to  other  gases,  depends  on  the  total  pressure  as  well  as 
the  amount  of  ozone  present.  When  the  absorption  due  to  ozone  takes 
place  at  a  uniform  pressure,  the  absorption  is  simply  a  function  of  the 
amount  of  ozone  m,  and  the  uniform  pressure  p,  along  the  path  (the 
temperature  effect  is  small  and  will  be  neglected  here).  Laboratory 
measurements  of  the  9.6fJ.  ozone  band  absorption  along  uniform  pres  - 
sure  paths  were  made  by  Summerfield  (1941)  and  by  Walshaw  (1957). 
Walshaw  expressed  the  absorption  in  terms  of  integrated  absorptivity 
over  the  entire  band  which  is  called  by  him  the  band  area,  and  the  ef¬ 
fective  band  width  by  some  others;  the  unit  being  wave  number  per  cm. 
Walshaw' 3  laboratory  measurements  will  be  used  as  the  basis  of  the 
following  analysis,  and  his  empirical  result  is  specifically  written 

A  =  W(m,  p)  (1) 

where  A  is  the  integrated  absorptivity  along  a  uniform  pressure  path, 
and  W(m,  p)  is  the  empirically  determined  function  shown  in  Fig.  1. 

In  the  case  of  absorption  of  solar  radiation  in  the  same  band 
by  atmospheric  ozone,  the  absorption  takes  place  in  a  path  along  which 
atmospheric  pressure  is  not  uniform,  and  the  resulting  absorption  de¬ 
pends  on  the  ozone  distribution  with  respect  to  pressure. 

There  are  several  ways  of  representing  the  vertical  distribution 
of  ozone.  The  commonly  used  ozone  density  (ozone  mass  per  unit  geo¬ 
metrical  volume)  is  inconvenient  for  theoretical  manipulation.  The 
ozone  mixing  ratio,  which  varies  by  many  orders  of  magnitude  through 
the  ozone  layer,  is  not  particularly  suitable  for  a  graphical  representation 
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or  numerical  calculation.  We  shall  use,  as  the  primary  quantity  de¬ 
scribing  the  vertical  ozone  distribution,  a  pseudo  ozone  density  defined 
on  a  logarithmic  pressure  scale.  The  exact  definition  is  given  in  the 
following. 

In  order  to  retain  the  advantage  of  using  atmospheric  pressure 
as  the  vertical  coordinate,  we  introduce  a  pressure -height  h,  defined  by 

h  =  -  Hte{p/p^)  .  (2) 

where  p  is  the  atmospheric  pressure,  p^  is  the  pressure  at  the  ground, 
and  H  is  a  constant,  having  the  dimension  of  length.  Although  the  constant 
H  can  be  chosen  completely  arbitrarily,  it  is  assumed  here  that  H  is  the 
scale  height  corresponding  to  a  temperature  of  220*K.  Then  h  is  only 
approximately  equal  to  geometric  height.  However,  since  the  difference 
between  the  pressure -height  and  geometric  height  is  not  of  much  im¬ 
portance  here,  we  may  simply  call  h  "height".  The  pressure  p,  now 
is  considered  to  be  a  function  of  the  new  independent  variable  h,  i.  e.  , 

p(h)  =  p^  exp(-h/H) .  (3) 

We  define  the  cumulative  ozone  distribution  function  m{h),  as 
the  amount  of  ozone  in  a  vertical  column  of  unit  cross  section  between 
the  ground  and  an  arbitrary  level  at  h,  and  also  define  the  pseudo  ozone 
density  u(h),  by 

u(h)  =dm/dh.  (4) 

If  m(h)  is  measured  in  cm  STP  and  h  in  km,  the  unit  of  u(h)  is 
cm  STP  per  km.  Because  of  the  above  choice  of  H,  the  pseudo  density, 
u(h),  is  nearly  equal  to  the  true  ozone  density  in  the  lower  stratosphere. 
The  true  density  at  any  height  is  given  by  (H/H)u(h)  ,  where  H  is  the 


scale  height  corresponding  to  the  temperature  at  the  height  under  consid¬ 
eration.  The  ozone  mixing  ratio  can  easily  be  shown  to  be  proportional 
to  u(h)/p(h). 

When  any  particular  ozone  distribution  is  given,  the  integrated 
absorptivity  of  the  9.6|J.  band  due  to  that  part  of  the  ozone  along  a  path 
between  the  ground  and  any  level  h,  which  is  denoted  by  A(h),  is  an  in¬ 
creasing  function  of  h.  However,  A(h)  is  a  function  of  the  ozone  dis¬ 
tribution  itself,  in  the  sense  that  A(h)  depends  on  all  values  of  u{h*)  in 
the  interval  0  <  h'  <  h.  No  laboratory  measurements  of  the  absorption 
along  such  a  variable  pressure  path  have  been  made,  and  theoretical 
calculations  of  A(h)  based  on  the  detailed  spectrum  of  the  9. bp  band 
are  not  practical.  However,  it  was  shown  by  Curtis  (1952)  that  in 
the  case  of  water  vapor  in  the  atmosphere,  the  absorption  along  a  vari¬ 
able  pressure  path  could  be  approximated  by  the  absorption  along  a 
uniform  pressure  path  of  the  same  amount  of  absorber  at  the  mean 
pressure  of  the  mass  distribution  of  the  absorber.  A  similar  method 
of  approximation,  including  also  the  temperature  effect,  was  independ¬ 
ently  developed  by  Godson  (1955).  Although  Curtis  indicated  that  the 
approximation  is  less  accurate  in  the  case  of  atmospheric  ozone  whose 
typical  vertical  distribution  is  quite  different  from  that  of  water  vapor, 
the  use  of  the  Curtis -Godson  approximation  was  recommended  by  the 
International  Ozone  Commission  (1956)  on  the  ground  that  no  better  ap¬ 
proximation  was  available.  Under  this  approximation,  we  have  for  a 
vertical  path 

A(h)  =*W[m(h),q(h)J  (5) 


Here,  the  two  arguments  of  Walshaw's  empirically  determined  function 
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W(m,p)  have  been  replaced,  respectively,  by  the  two  functions  of  h: 
the  amount  of  ozone  in  the  layer  below  the  level  h  which  is  given  by 

h 

m(h)  =  /  udh  ,  (6) 

o 

and  the  mean  pressure  of  the  ozone  distribution  in  the  same  layer  which 
is  defined  by 


q(h)  = 


h 

/  pu  dh 
o 


(7) 


If  an  exact  equality  is  to  be  assumed,  (5)  should  be  considered 

as  the  definition  of  the  spectroscopic  mean  pressure,  q  (h),  namely 

s 

A(h)  =  W[m(h).  qg(h)]  (8) 

The  Curtis -Godson  approximation,  as  introduced  by  the  authors,  meant 
that 

W(m(h),  qg(h)]  «  [m(h),  q{h)]  (9) 

However,  in  the  infrared  method  as  recommended  by  the  International 
Ozone  Commission,  a  tacit  assumption  is  made  that 


qg(h)“q(h)  (10) 

is  a  good  approximation  in  addition  to  (9).  The  approximation  of  this 
form  will  be  called  the  inverted  Curtis -Godson  approximation,  and  its 
accuracy  will  be  discussed  in  Section  4.  Meanwhile,  it  may  be  stated  here 
that  that  validity  of  approximation  (10)  does  not  necessarily  follow  from 
that  of  (9)  and  that  the  former  is  much  poorer  than  the  latter  particularly 
when  the  whole  ozone  layer  (h  =  ®®)  is  concerned. 

Atmospheric  ozone  at  every  height  emits  radiation  in  the  9.6^ 
band  according  to  its  own  temperatui’e.  The  emission  intensity  in  the 


9.6ij.  band,  E(h),  which  is  due  to  ozone  in  the  layer  below  the  level  h, 
and  which  reaches  the  ground  vertically,  may  be  given  by 

E(h)  =  /B(h)  ^^dh  (11) 

where  B(h)  is  the  average  black  body  emission  intensity  over  the  9.6|x 
band  at  the  temperature  at  each  height  h.  The  vertical  distribution  of 
temperature  is  considered  to  be  known  either  from  radiosonde  obser¬ 
vation  or  by  assumption,  so  that  the  black  body  intensity  in  each  case 
is  a  known  function  of  h. 

In  the  case  of  oblique  incidence,  as  observations  are  generally 
made,  u(h)  and  m{h)  in  the  above  formulas  have  only  to  be  multiplied 
by  the  secant  of  the  zenith  angle  of  the  path.  For  the  sake  of  simplicity, 
but  without  loss  of  general  validity  of  the  results,  it  is  assumed  that  ob¬ 
servations  are  made  with  vertical  incidence.  The  three  observational 
data  from  which  a  vertical  ozone  distribution  is  to  be  determined  are 
the  total  amount  of  ozone  m{<»),  the  integrated  absorptivity  of  the  9.6p. 
ozone  band  A(oo),  and  the  emission  intensity  E(“)  due  to  ozone  in  the 
entire  column  from  the  ground  to  the  top  of  the  atmosphere.  The  spectro¬ 
scopic  mean  pressure  of  the  entire  column,  q  (»),  is  immediately  deter- 

s 

mined  from  the  observed  values  of  m(oo)  and  A(oo)  by  definition  (8). 

The  problem  is  formulated  so  as  to  determine  the  function  u(h), 
from  the  three  integrals; 

The  total  ozone 

oo 

/u{h)dh  =  m{oo)  (12) 

o 

the  mean  pressure  under  approximation  (10) 
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oO 

/p(h)  u(h)  dh  oem{oo)  q  (oo)  (13) 

o  ® 

and  the  emission  intensity 

jB(h)  ^^dh  =  E(«)  (14) 

where  A(h)  is  related  to  u(h)  by  (5),  (6),  and  (7).  It  is  evidently  im¬ 
possible  to  determine  a  unique  solution  u(h)  from  a  finite  number  of 
integrals.  What  we  may  attempt  is  to  restrict  ourselves  to  seek  a  solu¬ 
tion  within  a  preselected  class  of  functions  which  contain  only  three 
independent  parameters,  and  to  solve  for  the  parameters. 

According  to  the  recommendation  of  the  International  Ozone 
Commission,  the  atmosphere  is  divided  into  four  layers.  The  ozone 
distribution  in  the  top  layer  above  the  10  mb  level  is  assumed  known. 

The  ozone  mixing  ratio  (consequently,  u(h) )  in  each  of  the  three  lower 
layers  is  assumed  to  be  proportional  to  a  given  function  of  pressure. 

The  parameters  of  the  ozone  distribution  in  this  model  are  the  amounts 
of  ozone  in  the  three  lower  layers.  The  assumed  linear  dependency  of 
u(h)  on  those  parameters  facilitates  calculation  of  a  solution.  The  re¬ 
sulting  ozone  density  profile,  however,  generally  shows  discontinuities 
at  the  boundaries  of  the  layers,  so  that  some  degree  of  subjectivity  is 
left  when  the  ozone  density  curves  are  to  be  smoothed.  Vigroux  (1959) 
modified  the  above  scheme  to  obtain  a  continuous  distribution  objectively. 
The  number  of  parameters  in  Vigroux' s  model  is  six,  although  only  three 
of  these  are  independent,  and  the  numerical  procedure  used  in  solving 
the  equations  is  rather  com.plex.  There  may  be  other  models,  but  it  is 
not  our  purpose  to  discuss  which  model  is  best.  Choice  of  a  model,  i.  e.  , 
selection  of  the  type  of  functions  as  possible  solutions  must  be  made 


inevitably  by  a  subjective  judgment  (or,  by  an  objective  judgment  based 
on  subjectively  chosen  rules)  on  how  well  the  model  fits  the  whole  var¬ 
iety  of  real  distributions  of  atmospheric  ozone.  Whichever  model  is 
chosen,  however,  it  is  not  possible  to  reproduce  most  details  of  the 
vertical  ozone  distribution.  Only  a  broad  profile  may  be  obtained  by 
adjusting  the  three  independent  parameters  in  the  model.  One  might 
expect  that  this  limit  of  resolution  would  be  improved  by  taking  infra¬ 
red  observations  at  different  zenith  angles  and  thus  by  supplying  more 
than  three  integrals  involving  the  unknown  function  u(h).  However, 
since  the  accuracy  of  the  infrared  observations  must,  as  is  shown 
in  the  following  section,  be  about  one  order  of  magnitude  better  than 
the  desired  accuracy  of  the  determined  ozone  distributions  in  the  strato¬ 
sphere,  observations  taken  at  different  zenith  angles  can  hardly  be  con¬ 
sidered  as  independent  information,  unless  a  marked  improvement  is 
introduced  in  correcting  the  raw  observational  data  for  the  influence  by 
absorbers  other  than  ozone. 

Therefore,  if  the  low  limit  of  vertical  resolution  obtainable 
by  the  infrared  method  is  accepted,  any  reasonable  model  may  be  se¬ 
lected.  Then,  when  the  observational  data  are  supplied,  the  vertical 
ozone  distribution  can  be  "uniquely"  determined. 

A  serious  question  still  remains  concerning  the  accuracy  of  an 
ozone  distribution  thus  determined,  if  inevitable  errors  in  the  observa¬ 


tional  data  are  taken  into  account. 
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3.  Error  analysis 

3. 1  Mathematical  development 

Assume  that  u(h)  represents  the  true  vertical  ozone  dis¬ 
tribution  and  that  m.(oo),  A(oo)  ,  and  E(oo)  are  the  basic  observational 
data  which  would  be  obtained  for  the  ozone  distribution  u(h)  if  there 
were  no  observational  errors.  Suppose  that,  for  some  reason,  we  ob¬ 
tain  from  actual  observations  slightly  different  values  m(oo)  +  6m(oo), 
A(oo)  +  6A(oo),  and  E{«>)  +  6E(<o),  which  would  be  the  correct  data  if  the 
ozone  distribution  were  u(h)  +  6u(h).  We  will  assume  that  the  relative 
errors  in  the  observed  data  are  sufficiently  small  so  that 
|6u(h)/u(h)|  <<  1  for  all  h. 

The  equation  for  total  ozone  is  now  written 
00 

/  (u  +  6u)dh  =  m(oo)  +  6m(oo) 
o 

or,  because  of  (12), 

^10  00 

/  6udh  =  6m{oo)  -  J  6udh  (15) 

°  ^10 

where  h^^  is  the  height  of  the  10  mb  level.  As  described  in  the  pre¬ 
ceding  section,  the  ozone  distribution  above  the  10  mb  level  is,  in  the 
infrared  method,  assumed  to  be  known.  If  the  assumed  distribution  is 
not  correct,  the  integral  on  the  right-hand  side  generally  does  not 
vanish,  and  this  error  in  the  assumption,  along  with  the  observational 
error  6m(oo)  ,  is  compounded  into  an  erroneous  estimate  of  the  ozone 
distribution  in  the  layer  below  the  10  mb  level.  Although  we  recognize 
the  different  causes  of  these  terms  on  the  right-hand  side,  they  are 
practically  inseparable.  Therefore,  we  may  simply  absorb  the  second 
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term  into  the  first  by  writing 

^10 

/  6udh  =  6m(oo)  (16) 

o 

From  the  definition  given  by  (7),  we  have 

CO 

/pudh  =  m(<io)  q{«>) 
o 

However,  q(<»)  is  not  obtainable  from  observations.  Instead,  q  (<»)  is 

s 

determined  from  the  observations  and  substituted  for  q(°°)  under  the  in¬ 
verted  Curtis -Godson  approximation  (10).  Therefore,  even  if  no  obser¬ 
vational  error  exists,  a  certain  error  is  made  in  estimating  q('=o).  In¬ 
stead  of  the  approximate  equality  (10),  we  may  specifically  write 

+  6j,q(‘")  (I'i') 

where  6^q(oo)  is  the  error  due  to  the  use  of  approximation  (10).  Since 

q  («o)  is  obtained  from  A(  oo)  and  by  the  definition  given  in  (8), 

s 

any  error  in  the  latter  quantities  will  give  rise  to  an  error  in  qg(oo);i.  e.  , 

A(oo)  +  6A(oo)  =  W[m(<»)  +  6m(oo),  q  (oo)  +  6q  (»)] . 

s  s 

Expanding  the  right-hand  side  in  a  Taylor  series  and  neglecting  higher 
order  terms  of  small  quantities,  we  obtain 

■  (Si'”'”’] 

where  the  subscript  oo  indicates  that  the  partial  derivatives  are  evaluated 
at  m  =  m(oo)  and  q  -  q(oo-).  (Since  we  consider  only  the  first  order  terms, 
the  difference  in  the  values  of  the  derivatives  at  q  =  q  (oo)  and  q(oo)  is 
neglected. )  Combining  both  types  of  error,  we  now  have  in  lieu  of  (13) 
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CO 

/  p(u  +  6u)dh  =  [m(oo)  +  5m('=o)]  (q(«>)  +  6  q(‘»)  +  6q  (“o)] 
o  c  s 

or,  again  neglecting  higher  order  terms  and  substituting  (18),  into  the 
expression  above  we  have 
^‘lO 

/  p  6u  dh  =  [q(oo  )  R(<»  )/S(oo )  6m{oo)  +  [1/S(oo)]  6A(co)  +  m(oo )  6  q(co)  (19) 
o 

where  R(<>o)  and  S(oo)  are  the  particular  values  at  h  =  oo  of  more  gene¬ 
rally  defined  functions 


R(h)  =  RMh),<,(h)l=(™  - 

and  (20) 

S(h)  .S[m(h),q(h)J.(^a 


These  derived  functions  are  shown  in  Fig.  2  as  functions  of  m  and 

p.  In  order  to  obtain  R  or  S  as  a  function  of  h,  it  is  necessary  only 

to  draw  a  curve  of  m  =  m(h)  and  p  a  q(h)  with  h  as  a  parameter. 

In  (19)  we  have  neglected  the  term 

oo 

/  p  6u  dh 


h 
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which  arises  for  the  same  reason  as  the  last  term  in  (15)  but  is  found 
by  conservative  estimate  to  be  minor  compared  to  the  right -ha.nd 
terms  of  (19). 

When  the  error  6E(oo),  in  the  observation  of  emission 
intensity  is  considered,  the  equation  corresponding  to  (14)  is  written 

After  integration  by  parts  and  rearrangemeiit  of  the  terms,  it  may 


be  written 


400 


Figure  2.  Derived  functions  of  W(m,  p)  :  (a)  R(m.p).  and  (b)  S(m,  p). 
For  definitions,  see  Eq  (20)  in  the  text. 


The  third  term  on  the  right-hand  side  depends  on  the  possible  error 
in  the  assumed  ozone  distribution  above  the  iO  mb  level,  and  also,  to 
some  degree,  on  the  error  in  the  computed  ozone  distribution  below  this 
level.  An  estimate  of  this  term  is  difficult  to  obtain  in  actual  cases, 
although  because  of  the  increasing  temperature  upward,  it  might  amount 
to  a  few  percent  of  E(oo).  We  shall  again  assume,  as  we  did  in  connection 
with  (15),  that  the  third  term  is  absorbed  into  the  first  term,  6E(<»). 

Since  the  Curtis -Godson  approximation  has  better  accuracy  in  esti¬ 
mating  the  absorption  than  in  estimating  the  mean  pressure,  the  approxi¬ 
mation  in  the  form  of  (9)  may  be  used  in  order  to  express  6A  of  the 
integrand  on  the  left-hand  side  of  (21)  in  terms  of  6u.  By  again  ne¬ 
glecting  higher  order  terms  in  the  Taylor  expansion,  we  have 

6A(h)  =  (^|^')  6m(h)+('l^')  6q(h)=R(h)  /  6udh+[S(h)/q(h>]  /  p6udh 

Integration  of  (21)  by  parts  after  substitution  of  6A  gives 

ho 

/  G(h)  5udh  =  6E(oo)  -  B(h,  )  6A(oo)  (22) 

o 

where  the  weighting  factor 

G(h)  .  -  R(h)  §  dh  -  p(h)  ]”  Ifi}  i  dl.  (23) 

can  be  numerically  computed  once  the  ozone  distribution,  u(h),  and 
temperature  distribution  are  specified. 

(16),  (19),  and  (22)  relate  three  different  weighted  integrals 
of  5u(h)  with  the  observational  and  theoretical  errors  on  the  right-hand 
sides.  It  is  not  possible  to  solve  these  equations  for  6u(h)  as  an 


arbitrary  function  of  h,  but  as  explained  before,  the  values  of  three 
parameters,  which  specify  a  solution  6u(h)  under  a  pre -selected 
model,  can  be  determined. 

The  model  which  is  chosen  here  for  numerical  examples 
(except  case  A')  is  similar  to  that  recommended  by  the  International 
Ozone  Commission.  The  atmosphere  below  the  10  mb  level  is  divided 
into  three  layers,  designated  by  Aj,  A^,  and  A^,  and  the  unknown 
parameters  are  the  errors  of  the  amount  of  ozone  in  each  layer: 

6x.  =  /  6udh  (i  =  1, 2,  3)  (24) 

^  A. 

1 

The  "true"  amount  of  ozone  in  each  layer 

X.  =  /  u  dh  (i  =  1 , 2,  3) 

^  A. 

1 

is  a  known  quantity,  since  u(h}  is  considered  here  to  be  given. 

As  to  variation  of  6u  with  respect  to  h,  6u  is  assumed 
within  each  layer  to  be  proportional  to  a  certain  power  of  pressure  p; 
namely,  in  each  A^, 

a.  /  a, 

6u(h)  =5  6x.  p  (h)//  p  dh  (25) 

^  /a. 

1 

A  description  of  the  layers  and  the  values  adopted  for  is  given  in 
Table  1  . 


Table  1  .  Details  of  the  model 

Layer,  i  Interval  A^  a. 

1  ground  (1000  mb)  to  tropopause  0 

2  tropopause  to  50  mb  -1 

3  50  mb  to  10  mb  0.8 
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Substitution  of  (25)  into  (16),  (19),  and  (22)  leads  to  linear  algebraic 
equations  for  6x. ,  which  may  be  more  conveniently  written  in  a  non- 
dimensional  form  by  considering  relative  errors,  6x7x^,  6m(oo)/m(oo), 
etc.  We  put  the  relative  error  of  the  amount  of  ozone  in  each  layer 

y.  a  6x./x^  (i  =  l,2,3),  (26) 

the  relative  observational  errors 


a  6m(oo)/m(oo),  z^  =  6A(«>)/ A(oo),  and  a  6E(oo)/E(«>),  (27a) 

and  the  relative  error  of  the  mean  pressure  due  to  the  use  of  the  in¬ 
verted  Curtis -Godson  approximation 


Z4  a  6^q(co)/q(oo)  (27b) 

We  may  call  Zj  the  input  error  of  type  j  (j=l,2,3,4),  and  y^^  the 
output  error  in  the  ith  layer  (i  a  1,2,  3). 

(16),  (19),  and  (22)  are  then  written  in  the  form 
3  4 

2^kiyi=I\j^j  (k  =  1.2,3)  (28) 

i=l  j=l 

wh  e  r  e 

a,  .  ax.  /m(oo) 
li  1  '  ' 


m(‘»)  q(oo)  /  p  dh 

A. 


1+a.  / 

.  ax.  /p  ^dh/ 

i 

a.  /r  a. 

.  a  X.  /  Gp  ^dh/  E(«>)  /  p  ^dh 

^  ^  A.  /  L  A.  J 


(29) 


for  i  =  1,2,3 


and 
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^1 1  "  ^24  ~ ^33  “  ^ 

*^12  "  *^13  "  ’^14  “  *^23  ~  ^ 
^21  ~  R(°°)/S(oo) 

^22  ~  A(M)/[S{«>)m(°o)q(oo)] 
=  -A(«>)B(h^Q)/E(«o) 


31 


(30) 


The  set  of  equations  (28)  can  be  solved  for  in  terms  of  z^;  i.  e.  , 

4 

yi  =  ^c.jZj  (i  =  1,2,3)  (31) 

j=l 

where  the  coefficient  represents  tbiC  sensitivity  of  the  error  of 
the  ozone  amount  in  the  ith  layer  to  the  input  error  of  the  type  j. 

In  order  to  examine  the  possible  change  in  the  sensitivity 
factors  due  to  different  modeling  of  5u,  we  consider  in  case  A',  be¬ 
low,  that  5u(h)  is  represented  by  four  connected  straight-line  seg¬ 
ments.  Assuming  that  6u(h)  is  constant  within  the  troposphere, 
zero  at  the  10  mb  level,  and  arbitrary  at  the  three  connecting  points- - 
tropopause,  80  ‘nb  and  50  mb,  we  may  reduce  (16),  (19),  and  (22) 
to  the  same  form  as  (28)  although  the  corresponding  definition  of 
coefficients  is  a  little  different  from  (29). 

3 . 2  Numerical  examples  and  discussion  of  the  result 

Since  the  sensitivity  factor  depends  on  the  temperature 
and  ozone  distributions  from  which  the  coefficients  in  (28)  are  calculated, 
various  combinations  of  these  distributions  should  be  examined  so  that 
a  general  conclusion  may  be  obtained  about  the  necessary  degree  of 
accuracy  of  the  observations  in  order  to  secure  a  desired  degree  of 
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the  determined  ozone  distributions.  However,  only  a  few  numerical 
examples  are  shown  below  in  which  the  assumed  temperature  and 
ozone  distributions  are  rather  typical.  Since  these  examples  disclose 
alarmingly  large  values  of  some  of  the  sensitivity  factors,  further  in¬ 
clusion  of  a  large  number  of  examples  is  not  essential  to  the  purpose 
of  this  study. 

The  assumed  distributions  of  ozone  and  temperature  are 
shown  in  Fig.  3  ,  and  their  combinations  for  each  case  are  given  in 

Table  2  .  Note  that  the  difference  between  A  and  A'  is  in  the  difference 
in  modeling  6u  as  discussed  above. 

Table  2  .  Combinations  of  ozone  and  temperature 
distributions  in  each  case  (see  Fig.  3 
for  the  distributions  which  are  referred 
to  by  numbers), 

Tremperature 


Case 

Ozone;  u(h) 

Troposphere 

Stratosphere 

A,  A' 

#1 

#1 

#1 

B 

#1 

#2 

#1 

C 

n 

U 

#1 

D 

#2 

#2 

#2 

The  final  results  of  the  sensitivity  factors  are  summarized 
in  Table  3.  Note  that  a  negative  value  of  c^^  implies  a  reduced  amount 
of  ozone  in  the  ith  layer  when  the  input  error  of  type  j  is  positive.  As 
an  example,  consider  Case  A  and  suppose  that  there  are  no  input  errors 
except  an  error  in  the  total  amount  of  ozone  (j  =  !)•  If  the  observed 
total  amount  of  ozone  is  larger  by  one  percent  than  the  true  value,  the 
amounts  of  ozone  in  the  troposphere  (i  =  1)  and  in  the  lower  strato¬ 
spheric  layer  below  the  50  mb  level  (i  =  2)  are  underestimated,  as  seen 
from  Table  3  ,  by  0.76  and  5.37  percent,  respectively.  On  the  other 


Figure  3.  Various  distributions  of  ozone  (a)  and  temperature  (b).  assumed  for  the  numerical 
examples.  See  Table  2  for  their  combination  in  each  case. 
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hand,  the  ozone  annount  in  the  upper  layer  between  50  and  10  mb 
(i  =  3)  is  overestimated  by  4.41  percent.  When  the  different  types 
of  input  errors  exist  simultaneously,  their  contributions  are  alge¬ 
braically  additive. 


Table  3  .  Sensitivity  factors,  in  various  cases, 

Out- 


put 

Case 

j  =  1 

2 

3 

4 

i 

error 

6m(oo) 

6A{“) 

6E(oo) 

m(oo) 

A(oo) 

E(oo) 

q(oo) 

A 

-0.76 

-  0.52 

3.25 

0,86 

6xi 

A' 

-0.64 

-  0.69 

3,  23 

0.83 

1 

B 

-1.54 

-  0.14 

4.76 

1.60 

^1 

C 

-0.69 

-  0.74 

3,62 

0.89 

D 

0.00 

-  3.00 

4.02 

0.00 

A 

-5.37 

18,27 

-4,82 

2.75 

5x- 

A' 

-4.41 

15.68 

-4.25 

2.33 

2 

B 

-4.17 

17.57 

-7.07 

1.62 

C 

-4.29 

16,08 

-4.63 

2.04 

D 

-5.17 

18.96 

-5.14 

3.17 

A 

4.41 

-  9.24 

2.13 

-1.48 

6x 

A' 

3.91 

-  7.91 

1.84 

-1.27 

3 

B 

3.88 

-  8.92 

3.13 

-0.98 

^3 

C 

6.18 

-14.80 

3.70 

-2.04 

D 

6.89 

-17.10 

4.11 

-2.94 

As  seen  in  the  rows  of  6x^/xj^  in  Table  3  ,  tropospheric  ozone 
is  distinctively  sensitive  to  the  error  of  emission  intensity.  There¬ 
fore,  if  the  amount  of  tropospheric  ozone  is  to  be  determined  within 
an  error  of,  say,  10  percent,  the  measurement  of  emission  intensity 
must  be  accurate  to  within  an  error  of  at  most  2  to  3  percent,  while 


the  other  types  of  input  error  may  be  larger  than  that.  This  finding 
agrees  in  general  with  Goody  and  Roach's  (1956)  contention  that  the 
infrared  method  (without  simultaneous  absorption  measurement)  is 
"capable  of  giving  useful  information  about  tropospheric  ozone". 
However,  there  is  the  possibility  as  in  Case  D,  that  the  accuracy 
of  the  absorption  measurement  is  almost  equally  important  in  deter¬ 
mination  of  the  tropospheric  ozone  amount,  although  the  assumed 
stratospheric  temperature  in  this  case  is  a  little  unusual. 

For  accuracy  of  the  determined  ozone  amount  in  the  second 
and  third  layers,  the  effect  of  the  emission  error  remains  at  about 
the  same  relative  level  as  for  the  tropospheric  ozone,  but  sensitivity 
to  the  other  types  of  input  error  is  increased  considerably.  In 
particular,  accurate  absorption  data  are  the  decisive  factors  in  ob¬ 
taining  reliable  ozone  distributions  in  the  stratosphere. 

If  the  ozone  amounts  in  the  two  stratospheric  layers  are  to 
be  determined  within  an  error  of  less  than  10  percent,  it  can  be  seen 
from  Table  3  that  the  error  of  the  absorption  data  must  be  less  than 
about  0.5  percent.  Although  the  above  error  analysis  has  been  made 
on  the  assumption  that  both  the  input  and  output  errors  are  sufficiently 
small,  the  result  may  at  least  qualitatively  be  extended  to  a  larger 
range  of  errors.  Then,  if  the  average  error  in  the  absorption  measure¬ 
ments  exceeds  5  percent,  the  stratospheric  ozone  distribution  deter¬ 
mined  from  such  observational  data  will  be  almost  totally  valueless. 
Furthermore,  since  the  absorption  error  affects  the  ozone  amounts 
in  the  second  and  third  layer  in  the  opposite  sense  from  each  other, 
even  a  slight  over  or  under  correction  of  the  raw  absorption  records 
for  atmospheric  water  vapor  could  result  in  deceptive  correlation 
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between  the  stratospheric  ozone  distribution  and  some  other  meteorolo¬ 
gical  parameters. 

Comparing  the  sensitivity  factors  for  Case  A'  with  those  for 
Case  A,  we  notice  that  their  magnitudes  are  reduced,  but  only  slightly, 
with  use  of  the  continuous  model  of  6u(h).  It  does  not  seem  likely  that 
a  substantial  reduction  of  sensitivity  can  be  achieved  by  proper  choice 
of  a  model,  although  some  models  may  prove  less  sensitive  than  others. 

Considering  the  infrared  method  as  a  means  of  determining 
stratospheric  ozone  distributions,  we  now  recognize  that  rigorous  care 
must  be  taken  to  reduce  the  magnitude  of  the  input  errors.  Since  there 
is  no  standardized  way  of  obtaining  the  absorption  and  emission  data  as 
in  the  case  of  the  total  amount  of  ozone,  assessment  of  errors  of  these 
observational  data  has  to  be  made  with  specific  instrumentation  and  with 
specific  procedures  through  which  raw  observational  records  are  trans¬ 
formed  into  numerical  values  and  also  corrected  for  influences  of  water 
vapor,  carbon  dioxide,  etc.  Discussion  of  these  technical  details,  how¬ 
ever,  is  beyond  the  scope  of  this  paper. 

An  estimate  of  the  fourth  type  error,  which  is  caused  by  the 
use  of  the  inverted  Curtis -Godson  approximation,  is  given  in  the  following 
section.  A  likely  value  of  6^q(<»)/q(«>)  is  found  to  be  about  minus  20  to 
30  percent,  although  it  depends  on  the  vertical  ozone  distribution  and  the 
air  mass  in  the  direction  of  observation.  If  the  values  of  sensitivity 
factors  in  Table  3  are  used,  then  it  is  found  that  this  error  alone  would 
rather  consistently  result  in  underestimation  of  the  ozone  amounts  in 
the  first  and  second  layers  by  about  20  and  50  percent,  respectively, 
and  in  overestimation  of  the  ozone  amount  in  the  third  layer  by  about  40 
percent.  Thus,  some  kind  of  correction  for  this  type  of  input  error  must 
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be  devised.  A  crude  way  of  eliminating  a  gross  systematic  error  in 
the  mean  pressure  is  to  increase  the  observed  spectroscopic  mean 
pressure  by  adding  a  constant  amount  or  multiplying  by  a  constant  fac¬ 
tor  for  all  observations.  A  flat  increase  of  1  0  mb,  which  Goody  and 
Roach  (1956)  adopted  for  practical  reasons,  may  be  justified  because 
such  a  procedure  should  have  reduced  a  large  portion  of  the  systema¬ 
tic  error  in  the  tropospheric  ozone  amount.  However,  a  more  refined 
method  of  correction  should  be  developed  for  determination  of  the  strato¬ 
spheric  ozone  distribution  because  of  the  higher  sensitivity  there, 

4.  On  the  inverted  Curtis -Godson  approximation 

4.1  The  case  of  a  single  Lorentz  line 

As  mentioned  previously  in  Section  2,  the  Curtis -Godson 
approximation  in  the  original  form  is  to  estimate  the  absorption  due 
to  an  absorbing  gas  which  is  distributed  with  different  pressure  along 
the  path  of  radiation.  This  is  achieved  by  considering  a  path  with  the 
■same  amount  of  the  absorbing  gas  at  a  uniform  pressure  which  is  equal 
to  the  mean  pressure  of  the  mass  distribution  of  the  absorbing  gas  on 
the  original  variable -pressure  path.  However,  in  the  infrared  method 
of  determining  the  vertical  distribution  of  the  atmospheric  ozone,  it  is 
necessary  to  estimate  the  mean  pressure  of  the  .ozone  distribution  from 
the  observed  absorption  along  a  variable -pressure  path.  A  simple 
method  of  doing  this  is  to  use  the  Curtis -Godson  approximation  in  an 
inverted  way;  that  is,  the  mean  pressure  is  taken  to  be  equal  to  the 
spectroscopic  mean  pressure,  the  latter  being  defined  as  the  pressure 
at  which  a  uniform -pressure  path  with  an  equal  amount  of  the  absorbing 
gas  would  give  exactly  the  same  absorption  as  observed  along  the  variable- 


pressure  path. 
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Unfortunately,  however,  the  accuracy  of  the  inverted  approxi¬ 
mation  is  not  as  good  as  one  might  expect  from  the  favorably  checked 
accuracy  of  the  original  approximation.  We  shall  first  demonstrate 
this  for  the  absorption  due  to  a  single  Lorentz  line.  The  extreme 
simplification  of  spectrum  has  to  be  made  for  mathematical  expediency. 
An  extension  of  the  result  to  the  absorption  due  to  the  9. bp  ozone  band 
will  be  shown  later  in  this  section. 

The  absorption  coefficient,  k(v,p),  of  a  Lorentz  line  is 

given  by 

k(v,p)=|^-2-^  (32) 

V  +  a 


where  S  is  the  line  intensity,  v  is  the  distance  from  the  line  center 
on  the  wave  number  scale,  and  a  is  the  half  width  of  the  line,  which 
is  assumed  as  usual  to  be  proportional  to  the  pressure,  p,  so  that  if 
is  the  half  width  at  a  standard  pressure  p^, 


a  =  ao(p/Po)  . 

Consider  a  distribution  of  absorbing  gas  with  respect  to  pressure, 
which  may  be  represented  by  the  cumulative  distribution  function  m(p) 
along  a  path.  The  integrated  absorptivity  (the  equivalent  width)  due  to 
the  entire  path  is  a  functional  of  the  distribiitioh  function  m(p)  and  is 
given  by  the  integral 


A[m(p)] 


I  {l  -  (' 

-00  ^  ^ 


M 

/  k(v,  p)dm)  ^  dv 

o  /  J 


(33) 


where  M  is  the  total  mass  (per  unit  cross  section)  of  absorbing  gas 
along  the  path  and  p  in  the  integrand  is  considered  implicitly  as  a 


function  of  rn. 
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where 


and 


In  the  assumed  case  of  a  single  Lorentz  line,  we  may  write 

OO 

Alm(p)]  =  2  /  {l  -  exp(-2-YN{v)  )}  dv  (34) 

o 


■y  =  SM/ (2it  a^) 


N(v)  =  / 


o  (v/a)  +1 


d(m/M) 


We  may  consider  a  special  case  in  which  the  whole  mass  of 
the  absorbing  gas  is  at  a  certain  pressure  p,  i.  e.  ,  the  case  of  uni¬ 
form  pressure  path.  The  integrated  absorptivity  of  such  a  path  is  a 
function  of  M  and  p,  and  it  is  denoted  specifically  by  A{M,  p). 

Integral  (34)  for  a  uniform  pressure  path  can  be  expressed  in  terms 
of  Bessel  functions  (after  Ladenburg  and  Reiche;  see  Elsasser,  1942), 
namely 

A(M,  p)  =  Zir  a^r  f('Y/r)  (35) 

where 

r  =  P/P^ 

and  the  function  f(x)  is  defined  by 


f(x)  =  xe  J^(ix)  -  ijj(ix) ) 

whose  numerical  values  are  listed  in  the  paper  by  Elsasser,  It  is 
easily  shown  that,  for  large  values  of  y, 

—  1  /Z 

A(M,  p)  ~  2a^(2Tr  V  r)  (36) 

and  for  small  values  of  y, 

A(M,  p)  «  2a^iT  \[1  -  y/{2r)  +  •  •  •  ] 


(37) 
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Integral  (34)  is  generally  impossible  to  solve  in  a  closed  form 
for  an  arbitrary  distribution  m(p).  The  only  exception  is  the  case  where 
rn(p)  is  a  linear  function  of  pj  that  is,  where  the  mixing  ratio  of  the  ab¬ 
sorbing  gas  to  the  air  is  constant  through  a  layer.  When  such  a  layer  has 
a  finite  thickness,  the  integral  (34)  may  be  expressed  in  terms  of  a  hyper - 
geometric  function  (Strong  and  Plass,  1950).  In  particular,  for  an  infinite 
layer  extending  from  the  ground  at  pressure  P^  1°  p  =  0,  the  cumulative 
distribution  function  with  a  constant  mixing  ratio  is  written 

na(p)  =  (pQ  -  P)l^/Po  =  (1  -  r)M  (38) 

and  the  integrated  absorptivity  is  given  by 

A[m{p)]  =  2a^Tr^r(y  +  i)  /^y)  (39) 

The  limiting  values  are,  for  large  values  of  y  • 

A[m{p)]  ~2a^(TrY)^  (40) 

and  for  small  values  of  y, 

A[m(p)  ]  =  2a^'tTY[  1  -  (2  In  2)y  +  ‘  •  *  ]  (41 ) 

The  spectroscopic  mean  pressure,  q  ,  as  defined  previously, 

s 

is  determined  from 

A(M,qg)  =  Afm(p)] 
or 

2a^TTr^f(Y/rJ  =  2a^-ir^r(Y  +  i)/r(Y)  (42) 

where 

r  =  q  /p 
s  ^s  ^o 


The  limiting  values  of  are  explicitly  determined  by  making  use  of  (36) 
and  (40),  or  (37)  and  (41);  i.  e.  , 
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q  “*0  /2  as  v  — ►  oo 

and 

'-►p^/(4  In  2)  =  0.361  as  y-*  °o  . 

The  mean  pressure,  q,  of  the  mass  distribution  (38)  in  this  ex¬ 
ample  is 

M 

q  =  /  pdm/M  =  p  /2  . 
o 

Thus,  the  inverted  Curtis -Godson  approximation  of  estimating  the  mean 
pressure  from  the  spectroscopic  mean  pressure  is  asymptotically  cor¬ 
rect  only  for  large  values  of  y.  At  the  other  extreme  end  of  y— ►O,  the 
spectroscopic  mean  pressure  is  about  30  percent  smaller  than  the  mean 
pressure.  The  ratio  q^/q  for  intermediate  values  of  y  is  shown  by 
the  top  curve  in  Figure  4  . 

It  may  be  shown,  on  the  other  hand,  that  the  original  Curtis - 
Godson  approximation  of  estimating  the  absorption  along  a  variable - 
pressure  path  is  exact  in  the  two  limiting  cases  of  y  =  0  and  y-^oo  , 
while,  for  intermediate  values  of  y,  the  approximation  overestimates 
the  absorption,  In  the  above  example  of  a  path  with  a  constant  mixing 
ratio,  however,  the  error  does  not  exceed  6  percent. 

As  the  second  example,  we  take  another  kind  of  simple  distri¬ 
bution  in  which  the  density  of  the  absorbing  gas  is  constant  throughout  a 
layer  between  the  ground  at  p^  =  1000  mb  and  at  a  top  level  at  pj^=10  mb. 
Since  the  effect  of  temperature  variation  with  height  on  the  air  density 
is  rather  sinall  in  the  atmosphere,  the  mixing  ratio  of  the  absorbing 
gas  may  be  assumed  to  be  inversely  proportional  to  the  air  pressure. 
Then  the  cumulative  distribution  function  is  given  by 
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m(p)  =  Mln(p^/p)/ln(p^/pP  (43) 

and  the  mean  pressure  of  this  distribution  for  the  entire  layer  is 

q  =  (Pq  -  Pi)/ln(p^/pj) 

=  215  mb 

The  spectroscopic  mean  pressure  is  determined  in  the  same 
way  as  in  the  first  example,  but  entirely  by  numerical  computation.  The 
result,  as  the  ratio  qg/q,  is  shown  by  the  bottom  curve  in  Figure  4 
The  ratio  in  the  limit  as  is  about  0.35. 

The  distributions  that  have  been  considered  in  the  above  ex¬ 
amples  do  not  correspond  to  any  typical  distribution  of  the  atmospheric 
ozone.  The  latter  has  a  principal  maximum  of  density  somewhere  be¬ 
tween  the  100  mb  and  10  mb  levels,  and  the  density  in  the  troposphere 
is  smaller  than  the  maximum,  by  about  one  order  of  magnitude.  There¬ 
fore,  a  major  part  of  the  atmospheric  ozone  concentrates  in  .a  smaller 
range  of  pressure  variation  than  that  given  by  the  constant  density  dis¬ 
tribution  in  the  second  example.  In  this  sense,  the  bottom  curve  in 
Figure  4  may  safely  be  considered  to  present  a  lower  bound  to  the 
ratio  qg/q  which  may  be  found  for  an  actual  ozone  distribution.  On  the 
other  hand,  the  top  curve,  corresponding  to  a  distribution  with  constant 
mixing  ratio  (i.  e.  ,  exponential  decrease  in  density),  may  present  an 
upper  bound. 

As  the  last  example  (still  assuming  a  single  Lorentz  line),  we 
consider  a  distribution  that  is  similar  to  that  of  atmospheric  ozone.  In 
order  to  simplify  the  computation,  a  continuous  distribution  with  respect 
to  pressure  is  replaced  by  the  superposition  of  four  uniform-pressure 
paths  at  different  pressures.  The  assumed  numerical  values  in  the 
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computation  are  shown  in  Table  4 

Table  4  .  Assumed  ozone  distribution 

Layer  Pressure  Amount  of  ozone  relative 
_ _  to  the  total  (M) _ 

1  503  mb  5.5% 

2  108  27.4 

3  26  54.8 

4  5.7  12.3 

The  mean  pressure  of  this  distribution  is 

q  =  7.3  mb  . 

The  spectroscopic  mean  pressure  is  shown  again  as  the  ratio 
in  Figure  4  . 

4,2  Application  to  the  case  of  the  9.6p  band 

It  is  too  laborious,  if  not  impossible,  to  compute  an  integral 
like  (3)  with  the  actual  spectrum  of  the  9. 611  band  of  ozone.  If  the  in¬ 
tensity  and  half  width  of  all  lines  in  the  band  are  identical  but  the  line 
centers  are  located  at  random,  a  .simple  statistical  model  of  band  ab¬ 
sorption  may  be  applied  to  the  computation.  It  is  shown  in  such  a  case 
(Kaplan,  1953;  Godson,  1955)  that  the  essential  problem  in  evaluating 
the  integral  is  reduced  to  that  of  the  case  of  a  single  isolated  line.  How¬ 
ever,  when  we  consider  the  9.6fx  ozone  band  as  a  whole,  such  a  simple 
statistical  model  seems  hardly  to  be  applicable.  Nevertheless,  if  we 
try  to  extend  the  above  result  to  the  ozone  problem,  it  is  necessary  to 
have  at  least  an  estimate  of  the  order  of  magnitude  of  y,  which  might 
be  characteristic  of  the  absorption  due  to  the  ozone  band,  because  the 
whole  consideration  in  this  section  would  become  unnecessary  if  such 
a  characteristic  value  of  y  were  found  to  be  very  large. 

The  line  intensity  distribution  in  the  9. bp  band,  taken  from 
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the  theoretical  computation  by  Filler  (1951)  agrees  well  with  the  loga¬ 
rithmic  ogive  which  is  considered  by  Godson  (1955)  in  his  statistical 
model.  The  empirical  curves  of  the  absorption  band  area  of  the  9.6|a, 
band  (Fig.  1  )  are  found  to  fit  fairly  well  in  a  range  of  the  ozone  amount 
less  than  0.4  cm  STP  to  the  theoretical  curves  due  to  the  Godson  model, 
if  values  of  the  two  parameters  in  the  model  are  properly  chosen  at  each 
pressure  for  the  curve  to  be  fitted.  From  the  two  values  of  the  para¬ 
meters  at  each  pressure,  we  may  determine  a  kind  of  mean  line  inten- 

4  12 

sity.  (Specifically,  it  is  an  empirical  value  for  (SS^/SS^*  )  /2it 

which  is  a  generalization  of  S/2ira^  in  the  simpler  statistical  model 
with  constant  line  intensity  and  half  width.  The  subscript  t  is  to  iden¬ 
tify  individual  lines  and  the  summation  is  to  be  taken  over  all  the  lines 
in  the  band.  )  The  values  of  S/2TTa^  which  are  thus  determined  to  repre¬ 
sent  the  9.6fj.  band  are  shown  in  Table  5 

Table  5.  S/2iTa  for  the  9. 6 u.  band 
o 

p  (mb)  20  100  4a0_ 

S/2Tra^  (cm  STP)"^  0.44  0.69  1.38 

The  dependency  on  the  pressure  at  which  curves  are  fitted,  and 
the  fact  that  a  good  fit  between  the  empirical  curves  and  the  theoretical 
is  not  obtainable  for  the  ozone  amount  greater  than  0.4  cm  STP  are  prob¬ 
ably  because  the  assumption  of  the  random  line  distribution  in  the  statis¬ 
tical  model  is  not  fulfilled  in  the  actual  ozone  band.  Thus,  some  ambi¬ 
guity  still  remains  in  just  what  value  of  S/Zira.^  should  be  taken  for  a 
variable -pressure  path,  although  the  order  of  magnitude  is  fairly  well 
established. 

In  the  last  example  of  numerical  computation  of  the  spectro- 
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scopic  mean  pressure,  the  mean  pressure  of  the  distribution  is  71.3 
mb.  We  may  take 

S/Ztra^  =  0.6  (cm  STP)"^ 

as  a  rough  but  reasonable  value,  in  order  to  make  correspondence  be¬ 
tween  the  single  line  model  and  the  9.6p,  band.  Then,  if  the  total  amount 
of  ozone  in  a  vertical  column  is  0.3  cm  STP,  the  spectroscopic  mean 
pressure  for  various  air  mass  ,  at  which  the  absorption  measurement 
is  made,  is  obtained  as  shown  in  Table  6  . 


Spectroscopic  mean  pressure  for  the 
ozone  distribution  given  in  Table 

2  3  5  «> 


Table 

Air  mass 
M  (cm  STP) 
Y 

q^/q 

q  (mb) 


J 1.5 

0.3  0.45 

.18  .24 

.65  .70 

46.3  49.9 


0.6  0.9 

.36  .54 

.73  .77 

52.0  54.9 


1.5  “o 
.90  “o 
.83  1.00 

59.2  71.3 


The  above  numerical  result  is  derived  only  for  a  particular 
case  of  the  ozone  distribution  under  crude  assumptions.  However,  it 
seems  to  be  evident  that  the  use  of  the  inverted  Curtis-Godson  approxi¬ 
mation  seriously  underestimates  the  mean  pressure,  possibly  by  about 
20  to  30  percent.  Implication  of  this  error  in  the  infrared  method  of 
determining  the  vertical  ozone  distribution  was  already  discussed  in 
preceding  sections. 


5.  Conclusion 

The  infrared  method  as  described  in  Section  2  can  provide  us 
with  a  broad  profile  of  the  vertical  distribution  of  atmospheric  ozone 
from  the  observations  with  optical  instruments  at  a  ground  station,  of 
the  total  ozone  amount,  the  integrated  absorptivity  of  the  9.6p  ozone 
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band  in  the  direct  solar  beam  and  the  ozone  emission  intensity  in  the 
same  band.  However,  the  determined  ozone  distribution  is  found  to 
be  very  sensitive  to  errors  in  the  observations.  If  the  ozone  amount 
in  each  layer  of  the  two  stratospheric  layers  is  to  be  accurate  to  an 
error  of  less  than  10  percent,  for  instance,  the  error  in  the  absorption 
measurement  {including  its  correction  for  effects  of  other  gases)  must 
be  less  than  about  0.5  percent.  It  is  also  found  that  the  use  of  the 
inverted  Curtis -Godson  approximation  leads  to  a  considerable  amount 
of  systematic  error  in  the  determined  ozone  distribution,  its  magni¬ 
tude  reaching,  for  instance,  50  percent  in  the  lower  stratospheric  layer. 

The  question  as  to  how  accurate  the  determined  ozone  distri¬ 
bution  should  be  depends,  of  course,  on  intended  use  of  the  results. 

For  the  purpose  of  quantitative  estimation  of  the  stratospheric  motion 
with  the  use  of  ozone  as  a  tracer,  an  error  of  more  than  5  percent  in 
the  determined  ozone  density  is  perhaps  not  acceptable.  On  the  other 
hand,  for  the  purpose  of  Inferring  a  large  scale  pattern  of  the  strato¬ 
spheric  circulation  or  studying  the  correlation  between  the  ozone  dis¬ 
tribution  and  tropospheric  meteorological  parameters  based  on  monthly 
or  seasonal  statistics,  a  purely  random  error  of  10  or  even  20  percent 
in  the  determination  of  the  individual  ozone  distribution  might  be  tolerated. 
Even  if  this  were  the  case,  however,  one  would  have  to  eliminate  scrupu¬ 
lously  any  systematic  error  in  the  observation  data  which  is  correlated 
to  atmospheric  conditions. 
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PART  II 


Absorption  Processes  in  the  Stratosphere  and  Mesosphere 

by 

J.  London  and  C.  Prabhakara 

1 .  Introduction 

It  has  become  increasingly  recognized  that  atmozpheric  ozone 
plays  a  dual,  and  in  both  cases  important,  role  in  the  dynamics  and 
energetics  of  the  upper  atmosphere.  Because  of  its  relatively  long  life 
below  30  km  (which  includes  the  region  of  maximum  concentration)  ozone 
participates  in  the  general  circulation  of  the  upper  troposphere  and  lower 
stratosphere  in  a  very  direct  sense.  As  a  result,  studies  of  the  ozone 
budget  are  closely  linked  to  an  understanding  of  circulation  processes 
in  this  region  of  the  atmosphere.  Also,  because  ozone  absorbs  solar 
ultraviolet  radiation  very  strongly,  it  contributes  to  the  energy  sources 
and  subsequent  temperature  structure  of  the  region  from  20-80  km. 

The  absorption  of  solar  radiation  in  the  stratosphere  and  meso¬ 
sphere  is  due  almost  completely  to  the  strong  absorption  by  ozone  in  the 

9 

spectral  region  \  -^  3200  A.  The  net  result  of  this  absorption  of  solar 
energy  is  first  to  produce  a  quasi -stable  distribution  of  ozone  and  then, 
through  ozone  absorption  in  the  strong  Hartley  bands,  to  produce  the 
temperature  maximum  at  the  stratopause  (about  50  km).  The  following 
discussion  describes  calculations  of  the  photochemical  equilibrium  dis¬ 
tribution  of  ozone  as  a  function  of  height,  latitude  and  season;  and  the 
subsequent  absorption  of  solar  radiation  in  the  stratosphere  and  meso¬ 
sphere  assuming  the  computed  ozone  distributions  obtain.  Some  of  the 
limitations  of  the  results  are  also  discussed. 
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2.  The  photochemical  distribution  of  ozone 

The  first  complete  theory  for  the  presence  and  distribution  of 
ozone  in  the  upper  atmosphere  was  given  by  Chapman  (1930).  Since  that 
time  extensions  to  both  theory  and  calculations  have  been  made  by  Wulf 
and  Deming  (1936),  Craig  (1950),  Diltsch  (1956),  and  others.  According  to 
the  theory,  the  presence  of  atmospheric  ozone  results  from  a  combination 
of  processes  involving  photodissociation  and  recombination.  Ozone  is 
formed  mainly  by  photodissociation  of  molecular  oxygen  and  then  re¬ 
combination  between  atomic  and  molecular  oxygen.  Ozone  is  then  de¬ 
stroyed  directly  by  photodissociation  and  by  recombination  with  atomic 
oxygen.  The  chief  photochemical  reactions  involving  the  formation  and 
destruction  of  ozone  are  assumed  to  be 


Og  +  hv  -*  O  +  O 

(1) 

0+0+M-*  Og+M 

(2) 

O+Og+M  —  Og+M 

(3) 

Og  +  hv  —  Og  +0 

(4) 

^3 

O  +  Og  20g 

(5) 

hv  represents  the  photo  energy  necessary  to  dissociate  molecular  oxy¬ 
gen  or  ozone!  ,  kg,  kg  are  the  reaction  rate  coefficients  describing 
processes  (2),  (3),  and  (5);  and  M  represents  any  third  body  in  the 
three -body  collison  described  in  Eqs  (2,  3).  The  presence  of  M  is 
necessary  in  order  to  preserve  the  total  energy  and  momentum  of  the 
colliding  particles  in  the  system. 

If  we  assume  that  equilibrium  among  these  five  reactions  obtains. 
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we  can  write  two  equations 


dwi 

w 


24  '^2  +  4  > 


I3  - 


kgn^r^  =  0 


(6) 


and 

^  =  -  fang  +  kgn^n^n  -  k^n^ng  =0  (7) 

where  n,  n^^  ,  ng,  ng  are  the  number  density  of  air  molecules,  atomic 
oxygen,  molecular  oxygen,  and  ozone  respectively;  and  £g,  fg  are  the 
number  of  quanta  of  solar  radiation  absorbed  per  molecule  per  unit  time. 


f2  =  /  (r2q^d\.=  /  or2q^^exp{-(rgN2  -  o-gNa)  dX. 

V  X 


(8a) 


^3  =  /  (Tg q;^  dX  =  /  (Tg expl-tTg  Ng)  d\ 


(8b) 


where  <r  ,  cr_  are  the  absorption  cross  sections  of  0„  and  O,;  q  . 

3  a  d  j  oo\ 

represents  the  number  of  quanta  cm“®3ec"^  per  unit  wavelength  in  the 
solar  beam  available  at  the  top  of  the  atmosphere;  Ng  and  Ng  represent 
the  total  number  of  oxygen  or  ozone  molecules  above  some  level  and  in 
the  direction  of  the  solar  beam. 

The  equilibrium  conditions  described  in  (6)  and  (7)  can  be 
solved  to  give 


n 


3  " 


*^2 


(fs/f3)kn-  ng(l  -  x)^ 
T^gkg)  +  n'^fl  x[ 


1 

3 


(9) 


where  k  =  kg/kg  and  x  =  (ng/kn-iig). 

Since  any  atmospheric  gas  can  act  as  the  third  body  in  reactions 
(2)  and  (3),  Ng  (nitrogen)  and  ©2  will  be  the  third  body  99  percent  of 
the  time  except  at  levels  above  80  km  where  the  number  of  oxygen  atoms 
becomes  qtiite  large.  The  relative  efficiencies  of  N  and  0_  acting  as 

8  3 
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third  bodies  in  the  collision  processes  described  by  Eqs  (2,  3)  are  not 
very  well  known,  and  we  will  assume  that  both  gases  are  equally  efficient. 
For  a  mixed,  dry  atmosphere,  the  number  of  nitrogen  molecules  is  equal 
to  3,76  and  kn.n^  =  kng(4.76}  n^. 

Thus 

(10) 


.25n® 

3 


(f2/fa)k(l  -  x)^ 

1  -  x) 


In  order  to  simplify  this  relationship  further,  it  is  convenient  to  evaluate 
the  order  of  magnitude  of  the  various  terms  and  to  make  suitable  approxi¬ 
mations.  The  orders  of  magnitude  given  in  the  table  belov/  represent 
average  values  of  the  parameters  in  Eq  (10).  The  references  for  these 
values  are  discussed  below. 

We  see  that  for  all  levels  x  ^  10“^  1  and  we  can  write 


=  2.25n 


(f2/f3)k 

^Ik^2^+  ng 


(11) 


Also,  at  levels  below  about  50  km,  the  first  term  in  the  denominator  is 
smaller  than  the  second  term  by  a  factor  of  at  least  100.  Therefore,  for 
levels  below  50  km,  we  have 

ng  =  2.25n2[f2/f3(k)ng]^  (12) 

This  equation  is  the  same  as  that  given  by  Diitsch  (1956). 

According  to  Eq  (12)  the  equilibrium  ozone  concentration  below 
50  km  depends  on  the  solar  ultraviolet  ratio  fg/fg  and  the  recombination 
ratio,  k.  Changes  in  ng  would  require  changes  in  the  relative  values  of 
fp  as  compared  to  fg.  The  role  of  the  recombination  rate  ratio  k  is  to 
stabilize  the  ozone  production  in  the  region  below  50  km.  If  we  regard 
Eqs  (1-5)  as  the  primary  reactions  for  determining  the  photochemical 
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equilibrium  distribution  of  ozone,  we  require  the  following  data  to  com¬ 
pute  ns  from  Eq  (11)  or  Eq  (12). 

(a)  q  .  the  spectral  distribution  (in  absolute  intensities  of 

ccA. 

the  solar  beam  at  the  top  of  the  atmosphere  for  wavelengths 
covering  the  interval  of  Og  and  absorption  in  the  ultra¬ 
violet  and  visible.  For  purposes  of  the  present  computation, 
this  includes  the  region  1700  A  to  7500  A. 

(b)  ffg,  (Ta  ■  ■  absorption  cross  sections  of  Og  and 
in  this  spectral  interval. 

(c)  kj,  k.g,  kg-  -  the  reaction  rate  coefficients  appropriate  to 
Equations  (2),  (3),  and  (6). 

(d)  The  distribution  of  temperature  and  pressure  and  the  re¬ 
sulting  distribution  of  in  the  atmosphere  (20-  80  km). 

2. 1  The  spectral  distribution  of  solar  radiation 

The  solar  radiation  curve  approximates  an  equivalent  black 

body  radiator  (of  about  6000  °K)  in  the  visible  and  near  infrared  part  of 

the  spectrum.  Because  of  the  strong  absorption  by  ozone  it  was  impossible 

in  earlier  years  to  obtain  solar  spectra  that  extended  below  about  2800  A. 

Rocket  observations  taken  since  1946,  however,  have  extended  the  observed 

0 

spectrum  down  to  2200  A  and  in  some  cases  through  sections  of  the  far 

O 

ultraviolet.  The  origin  of  the  solar  spectrum  in  the  interval  1500  A  - 
2000  X  is  not  clearly  understood  and  there  has  been  some  speculation  of 
the  character  of  the  solar  emission  curve  for  this  region.  The  observations 

9  O 

show  that  the  solar  spectral  curve  in  the  region  3000  A  to  2000  A  falls 
off  sharply  to  an  equivalent  black  body  curve  of  about  5000 °K. 

The  energy  spectrum  used  in  the  present  calculations  were  taken 
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from  Johnson  (1954),  and  Johnson  et  al  (1954)  who  used  the  available  ob¬ 
servations  to  reconstruct  a  solar  energy  curve  for  the  range  2200  X  to 

O 

8000  A.  We  then  extrapolated  these  values  to  fit  an  equivalent  black  body 
curve  of  5000  “K  at  2000  A  to  4500  °K  in  the  wavelength  region  1900  A 
to  I  700  A. 

After  the  calculations  for  n^  were  completed,  additional  rocket 
spectra  became  available  (Detwiler  et  al,  1961)  indicating  that  the  inten¬ 
sity  distribution  in  the  solar  ultraviolet  at  175oX  corresponds  to  an 

equivalent  black  body  radiating  temperature  of  about  4700 'K.  Since  the 

1. 

equilibrium  concentration  of  ozone  below  50  km  is  given  by  n2a(f2/f2)^  an 
increase  in  the  value  of  f2  (holding  constant)  will  increase  the  ozone 
concentration.  This  correction  should  apply  to  the  ozone  distribution  at 

9 

the  levels  40-50  km  where  radiation  in  the  interval  1900 -2000  A  is 
most  effective  as  will  be  shown  below. 

2.2  Absorption  cross  sections  for  Og  and  Os 

I.  Oxygen:  There  are  three  important  regions  in  the  ultra¬ 
violet  for  absorption  by  molecular  oxygen.  These  are 

a.  Schumann -Runge  continuum  (1250  A  -  1750  A).  This 
region  i.s  largely  responsible  for  the  dissociation 

of  oxygen  above  80  km. 

b.  Schumann -Runge  bands  (1751  A  -  2026  A).  These  bands 
are  absorbed  in  the  mesosphere  principally  in  the 

layer  50-80  km. 

0 

c.  Herzberg  continuum  (X.<  2424  A).  Although  absorption 
in  this  region  represents  the  weakest  of  the  three, 
absorption  in  the  spectral  interval  2000  A  -  2200  A 


is  responsible  for  the  major  production  of  ozone  in 
the  stratosphere  and  mesosphere.  Solar  radiation 

®  a 

at  2000  A- 2200  A  penetrates  deeply  into  the  strato¬ 
sphere  resulting  in  dissociation  of  molecular  oxygen 
at  a  level  where  the  chance  for  three  body  recombi¬ 
nation  to  form  ozone  is  quite  good. 

The  oxygen  absorption  cross  sections  used  in 
Eq  (8a)  for  the  interval  1700  A-  1950  A  were  taken 
from  Watanabe  (1958).  For  the  interval  2I00A  - 

A 

2424  A  it  was  assumed  that  the  oxygen  absorption 
was  pressure  dependent  in  the  manner  discussed 
by  G6tz  (1951).  The  absorption  cross  sections  were 
approximated  for  zero  pressure  from  the  values  given 
by  G8tz.  For  the  spectral  region  1950A-2100A 
the  values  were  extrapolated  between  the  data  as 
discussed  above. 

2.  Ozone;  The  important  regions  for  ozone  absorption  s.re 
the  weak  Chappuis  bands  in  the  visible  spectrum,  the 
Huggins  bands,  and  the  contiguous  Hartley  bands.  Of 
these  three,  Hartley  absorption  is  by  far  the  most  impor¬ 
tant.  Although  the  ozone  absorption  cross  section  for 
X  <  2020 X  increases  again,  this  increase  is  in  a  region 
of  decreasing  solar  energy  and  represents  only  a  minor 
contribution  to  the  absorption  process  in  the  stratosphere 
and  mesosphere.  The  absorption  cross  sections  for  ozone 
were  taken  from  Vigroux  (1953)  for  wavelengths  larger  than 
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2300  A  and  from  Inn  and  Tanaka  (1953)  and  Watanabe  (1958) 
for  the  region  below  2300  A.  These  values  represent  the 

results  of  careful  laboratory  experiments  and  are  probably 

the  most  reliable  available  at  the  present  time. 

The  distribution  of  solar  radiation  at  the  top  of 
the  atmosphere  (q  .)  and  absorption  cross  sections 
((r2>(T^)  used  in  the  present  calculations  are  given  in 
Figs.  1(a),  1(b).  Curves  (1)  and  (2)  of  Fig.  1(a)  refer 
to  the  absorption  cross  sections  of  ozone  and  molecular 
oxygen  (cm®molec“^).  Curve  (3)  refers  to  the  solar  energy 
arriving  at  the  top  of  the  atmosphere  on  a  surface  perpen¬ 
dicular  to  the  sun's  rays  (units  of  right  hand  coordinate 
are  quanta  cm“®sec"^  A"^).  Fig.  1(b)  shows  the  absorption 
cross  section  for  ozone  [curve  (1)]  and  the  solar  energy 
distribution  [curve  (2)]  on  a  compressed  wavelength  scale. 

2.3  The  rate  reaction  coefficients 

The  recombination  rate  coefficients  k^  and  ks  were  taken  from 
recent  re -analysis  by  Campbell  and  Nudelman  (I960)  of  the  data  of  Gliss- 
mann  and  Schumacher.  They  give  as  best  estimates  for  the  individual 
rate  coefficients  for  the  reaction  given  by  Eqs  (3,  5); 

kg'  =  1.87  X  10^^exp[149/T]  cm®/gm  mol^/sec 
kg  =  5.44  X  10^®exp{-2133/T]cm®/gm  mol/ sec 

where  the  third  body  in  the  reaction  for  k^'  was  O^.  Campbell  and 
Nudelman  assume  that  the  efficiency  of  kg'  is  reduced  by  one -third 
when  the  third  body  in  Eq  (3)  is  Og  (or  presumably  Ng).  Therefore, 
we  can  write  for  the  reaction  rate  ratio  for  the  atmosphere! 


Fig.  1(a)  The  spectral  distribution  of  absorption  cross  section 

(cm®  molec"^)  for  ozone  and  oxygen,  and  solar  ratiation 
(quanta  cm“®sec~^  A~^)  received  at  the  top  of  the  atmosphere 
on  a  surface  perpendicular  to  the  sun's  rays.  Spectral 
interval  1700^-  2424  X. 


Fig.  1(b)  The  same  for  ozone  and  solar  radiation.  Spectral  interval 
1700  A  -  7500  A. 
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1.9  X  1 0"®^  exp(2282/T)  cm^/moleciile. 


(13) 


The  estimate  of  k  thus  derived  can  be  compared  with  earlier  values 
as  given  by  Eucken  and  Patat  {1936}  for  the  range  of  temperatures 
observed  in  the  stratosphere  and  mesosphere.  It  is  seen  from  Fig. 

2  that  the  two  curves  are  quite  close  in  the  region  of  relatively 
low  temperatures  (~’200*K)  and  differ  only  by  a  factor  of  2  at 
temperatures  corresponding  to  the  stratopause  {~290*K).  As 
will  be  indicated  below,  however,  variations  in  k  are  important  only 
at  the  level  of  the  mesopause  where  the  self -regulating  mechanism 
on  ozone  formation  is  weakened.  It  is  ia  this  region  that  the  dif¬ 
ference  between  the  two  curves  shown  in  Fig.  2  is  quite  small. 

The  equilibrium  distribution  of  ozone  at  levels  above 
about  60  km  depends  on  the  values  for  the  rate  coefficients  used 
in  Eq  (11).  Since  both  k^  and  kg  describe  three  body  recombination 
processes,  the  ratio  of  their  values  has  been  taken  as  independent  of 
temperature.lt  is  clear  from  Eq  (11)  that  as  the  ratio  increases, 

the  equilibrium  value  for  ozone  will  decrease.  The  value  used  for 
the  rate  ratio  in  the  present  computations  was  kj^/kg  =1,5.  This 
value  is  somewhat  low  when  compared  to  the  ratio  of  the  individual 
rate  coefficients  given  by  Bates  and  Nicolet  (1950)  and  Campbell  and 
Nudelman  (I960)  as  shown  in  Table  1  above.  The  use  of  the  larger 
rate  ratio  would  not  affect  the  equilibrium  ozone  concentration 
below  65  km  but  would  redxxce  the  values  above  75  km  by  about 
60  percent. 


^2^  ^3  (cm  /molec) 
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200  217  238  263  294  333 

T  ®K 

Fig.  2.  Variation  of  the  recombination  rate  ratio  kg/ka  with  temperature 

after  Campbell  and  Nudelman  (C  and  N)  and  Eucken  and  Patat  (E  and  P). 


2.4  The  model  atmosphere 


Molecular  oxygen  was  assumed  to  act  as  a  reservoir  for  the 
formation  of  atomic  oxygen  and  ozone.  Consequently,  the  distribution  of 
molecular  oxygen  below  about  80  km  is  a  function  only  of  the  distri¬ 
bution  of  pressure  and  temperature.  If  the  concentration  of  ng  is  taken 
to  be  21  percent  of  that  of  total  atmospheric  particles,  we  have 

1  p 

""  TTTB  kT 
=  1.52  X  10^®  p/T 

where  p  is  given  in  mb. 

There  are  sufficient  data  available  to  provide  for  realistic 
models  of  the  atmosphere  up  to  about  80  km.  For  the  present  compu¬ 
tations  the  distribution  of  pressure  and  temperature,  for  winter  and 
summer,  was  taken  from  the  atmospheric  cross  sections  given  by  Mur- 
gatroyd  (1957),  These  values  were  extrapolated  to  polar  and  equatorial 
latitudes  where  necessary. 

3.  The  equilibrium  ozone  distribution 

The  photochemical  equilibrium  distribution  of  ozone  was  com¬ 
puted  from  Eq  (11)  by  a  numerical  procedure.  It  was  assumed  that  the 
ozone  concentration  above  90  km  is  zero  and  that  the  total  depth  for  mole¬ 
cular  oxygen  above  90  km  is  known.  Then,  since  the  temperature  and 
pressure  is  specified  at  all  altitudes  and  latitudes,  the  molecular  oxygen 
concentration  and  the  reaction  rate  coefficients  are  known  for  each  level. 
The  computation  was  by  iteration  starting  with  the  spectral  distribution 
of  solar  radiation  at  the  top  of  the  atmosphere  and  computing  downward 
from  90  km  in  2 -km  steps. 
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The  computed  distributions  depend  somewhat  on  the  temperature 
distribution  in  the  atmosphere  and,  at  some  levels,  depend  rather  strongly 
on  the  zenith  position  of  the  sun.  In  order  to  study  the  effect  of  this 
latter  relationship,  the  photochemical  distribution  was  computed  for 
the  temperature  distribution  of  0*  latitude  for  various  solar  zenith 
angles.  The  results  are  shown  in  Fig.  3.  The  curves  marked  ^  =  0*, 

57",  86.75°  correspond  to  noon  at  the  equator  during  the  equinoxes, 
mean  zenith  angle  at  the  equator  at  summer  (winter)  solstice,  and  mean 
zenith  angle  at  60°  at  winter  solstice  respectively. 

It  can  be  seen  from  the  variation  of  the  curves  in  Fig.  3  that 
there  is  marked  solar  control  of  the  equilibrium  distribution  below  about 
45  km.  All  curves  show  the  well  recognized  pattern  of  increasing  ozone 
to  a  maximum  and  then  a  sharp  decrease  above.  The  equilibrium  maxi¬ 
mum  results  from  the  dissociation  of  molecular  oxygen  in  that  part  of 
the  Herzberg  continuum  where  the  ozone  absorption  is  a  minimum 
(2000  -  2200  A).  As  a  result,  solar  energy  can  penetrate  downward  to 
about  25  km  for  a  vertical  sun  and  o  about  35  km  for  a  highly  oblique 
sun.  Also,  for  large  optical  paths  (high  zenith  angles)  the  solar  radiation 
is  strongly  depleted  and  the  equilibrium  ozone  value  at  the  level  of  the 
maximum  concentration  is  down  by  a  factor  of  about  30  from  the  peak  at 
vertical  incidence  (3  X  10^^  cm"®  as  compared  to  10^® cm"®). 

The  general  shape  of  the  ozone  distribution  curves  as  well  as 
the  decrease  in  ozone  concentration  and  the  increase  in  the  height  of 
maximum  concentration  with  increasing  zenith  angle  follow  the  pattern  of 
a  typical  "Chapman  layer"  first  derived  to  describe  the  formation  of 
layers  in  the  ionosphere  (Chapman,  1931). 


The  photochemical  equilibrium  distribution  of  ozone  for 
solar  zenith  angles.  (The  temperature  distribution  corr 
to  that  for  latitude  0*.  ) 
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It  can  be  seen  that  in  the  region  45-55  km,  the  ozone  distribution 
is  independent  of  zenith  angle.  Above  this  level,  however,  increased 
path  lengths  through  the  oxygen  atmosphere  again  result  in  lower  ozone 
concentrations.  The  slight  bulge  in  the  ozone  distribution  that  appears 
near  70  km  is,  in  part,  due  to  the  strong  oxygen  dissociation  that  takes 
place  at  the  head  of  the  Schumann -Runge  band  system. 

Equilibrium  ozone  distributions  were  computed  for  various  lati¬ 
tudes  for  winter  and  summer  solstices.  These  are  shown  in  Fig.  4(a) 
for  summer  latitudes  0°,  30“,  60“,  and  80';  and  Fig.  4(b)  for  winter 
latitudes  O',  30*,  and  60'.  Since  the  equilibrium  concentrations  have 
long  half-lives  below  45  km,  where  the  solar  zenith  angles  become  im¬ 
portant,  in  each  computation  a  mean  zenith  angle  was  assumed  correspond¬ 
ing  to  the  average  angle  between  sunrise  and  sunset  for  the  winter  (sum¬ 
mer)  solstice.  The  average  zenith  angle  used  for  each  latitude  is  indi¬ 
cated  in  the  legend  in  Figs.  4(a),  (b).  The  seasonal  and  latitudinal  vari¬ 
ations  are  quite  marked,  particularly  below  40  km  and  above  60  km.  The 
largest  seasonal  variations  are,  of  course,  at  high  latitudes.  It  should  be 
noticed  that  at  50  km  the  computed  ozone  concentration  is  reasonably 
constant  at  approximately  5  X  lO’'^  molec/cm'®. 

At  30  km  the  ozone  concentration  at  60'  latitude  decreases  by  a 
factor  of  about  50  from  summer  to  winter  (but  only  by  a  factor  of  10  if 
the  peak  concentrations  are  compared).  At  70  km,  the  decrease  in 
photochemical  ozone  concentration  from  summer  to  winter  at  high  lati¬ 
tudes  is  about  fourfold. 

Although  the  variations  in  the  vertical  distributions  shown  in 
Fig.  4(a),  (b)  are  mostly  due  to  the  effect  of  solar  zenith  angle,  the  atmo¬ 
spheric  temperature  distribution  could,  through  the  temperature  sensitive 
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recombination  rate  coefficients,  also  affect  this  distribution.  The  most 
important  is  the  rate  ratio  k  =  k^/ka,  which  decreases  over  two  orders 
of  magnitude  for  an  increase  of  100‘C  in  atmospheric  temperature. 

The  temperature  effect  on  the  computed  vertical  ozone  distri¬ 
bution  is  shown  in  Fig.  ,  5(a).  The  computations  were  made  for  a  tem¬ 
perature  distribution  corresponding  to  30*  summer  and  30“winter.  In 
order  to  isolate  the  temperature  effect,  the  average  zenith  angle  ^  =  45* 
was  used  for  both  computations.  The  assumed  temperature  structure  is 
shown  in  Fig.  5(b).  Despite  the  fact  that  there  is  a  larger  temperature 
variation  below  50  km  than  above,  the  effect  on  the  computed  ozone  dis¬ 
tributions  is  negligible  below  50  km  and  results  in  about  a  20  percent 
reduction  above.  This  is  because  at  levels  below  50  km,  the  production 
of  ozone  is  very  sensitive  to  the  available  solar  energy  in  the  region 

O 

of  important  overlapping  oxygen  and  ozone  absorption  u*  »  2000-2400  A). 
Thus,  initially  higher  temperature  results  in  smaller  recombination  to 
form  ozone  which  gives  less  ozone  and  therefore  deeper  penetration  of 
the  solar  beam.  This  results  in  more  energy  being  made  available  to 
produce  slightly  more  ozone  in  the  level  below.  There  is,  therefore,  an 
efficient  self -regulating  mechanism  that  produces  rather  stable  equi¬ 
librium  ozone  concentrations  despite  possible  temperature  changes  at  or 
just  below  the  stratopause  (50  km).  The  consequence  of  this  stable  re¬ 
action  for  direct  heating  of  the  upper  stratosphere  has  been  discussed  by 
Craig  and  Ohring  (1958),  who  show  that  this  effect  will  inhibit  large  tem¬ 
perature  changes  in  the  vicinity  of  the  stratopause. 

In  the  region  above  60  km,  the  situation  is  somewhat  different. 

As  before,  higher  initial  temperatures  result  in  lower  ozone  concen¬ 
trations.  At  these  heights,  however,  the  dissociation  of  molecular  o'scygen 
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Fig.  5(a)  The  photochemical  equilibrium  distribution  of  ozone 
corresponding  to  the  temperature  structures  for  30' 
winter  and  30*  summer.  (Calculated  for  ^  =  45“) 
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is  produced  by  absorption  at  wavelengths  where  the  depletion  of  solar 
radiation  by  the  total  ozone  column  is  not  very  important  (i.  e. 

X  <  1900  A).  Lower  ozone  concentrations  above  60  km,  therefore,  have 
relatively  little  effect  on  the  available  energy  at  these  levels.  Thus,  the 
influence  of  higher  temperatures  such  as  found  at  high  latitudes  in  the 
mesosphere  during  winter  as  compared  to  summer,  particularly  when 
combined  with  the  effect  of  increased  zenith  angles,  result  in  decreased 
equilibrium  ozone  amounts  in  the  upper  mesosphere.  (Compare,  for 
instance.  Figs.  4(a)  and  4(b)  .  )  In  the  absence  of  any  self-regulating 
mechanism,  some  anomalous  heating  of  the  air  in  the  mesosphere  could 
then  result  in  a  decrease  of  the  total  ozone  above,  say,  60  km  allowing 
for  a  greater  penetration  of  near  ultraviolet  radiation  to  the  region  of 
the  stratopause. 

The  vertical  ozone  distributions  shown  in  Figs.  4(a),  (b)  are  very 
similar  to  those  given  by  Diitsch  (1956),  with  two  notable  exceptions. 

The  level  of  maximum  ozone  is,  in  all  cases,  higher  in  the  present  cal¬ 
culations,  and  at  the  level  of  maximum  ozone  concentration  and  above, 
the  present  calculations  give  larger  ozone  amounts  than  those  computed 
by  Ddtsch.  Below  the  maximum  level,  DQtsch's  values  are  larger. 

Both  of  these  differences  stem  primarily  from  the  differences  in  model 
atmospheres  (temperature  distribution)  assumed  and  the  present  ab¬ 
sorption  cross  sections  used  for  the  computations. 

4.  The  absorption  of  solar  radiation  by  atmospheric  oxygen  and  ozone 

Monochromatic  solar  radiation  in  the  spectral  region  (1700  - 

O 

3000  A)  is  depleted  as  it  penetrates  through  the  mesosphere  chiefly  by 
absorption  by  molecular  oxygen  and  ozone.  In  that  case,  the  transmission 
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of  solar  radiation  down  to  any  atmospheric  level  is,  from  Eq  (8a,  b),  as 

y^{Z)  =  q^(Z)/q^(«>)  =  expI-(cr^2N2(Z)  +a-j^3N3(Z)]  (14) 

The  depth  of  penetration  of  solar  radiation  in  the  atmosphere  was  com- 

•  S 

puted  for  the  spectral  region  1720  A  to  3 100  A  for  an  atmospheric  transmission  of 
5  percent  (i.  e.  ,  assuming  that  95  percent  of  the  solar  energy  is  absorbed 
in  the  layer  above  the  level  indicated)  and  for  various  latitudes  and 
seasons.  was  calculated  from  the  equilibrium  distributions  as  dis¬ 

cussed  above.  The  results  are  shown  in  Fig.  6  for  two  of  the  extreme 
cases,  those  of  30*  summer  and  60*  winter.  The  portion  of  the  curves 
lying  above  85  km  is  quite  unreliable  since  it  depends  strongly  on  the 
total  number  of  oxygen  molecules  above  this  level.  This  is  the  region 
where  oxygen  dissociation  is  quite  strong,  and  there  is  probably  large 
variation  in  the  molecular  oxygen  concentration  due  to  vertical  motions 
and  associated  recombination  processes. 

The  difference  in  the  depth  of  penetration  shown  by  the  two 
curves  is  roughly  constant  at  about  10  km.  This  means  that,  on  the 
average,  the  seasonal  and  latitudinal  variations  of  the  penetration  of 
solar  radiation  can  be  as  large  as  about  10  km.  The  difference  in  pene¬ 
tration  through  the  strong  Schumann -Rung  e  bands  due  to  molecular 
oxygen  absorption  is  quite  evident.  It  is  seen  that  some  solar  energy 
can  filter  down  through  the  absorption  minima  to  levels  of  about  60  km. 

The  most  important  region  as  far  as  ozone  production  is  con- 
cerned  is  the  region  of  absorption  overlap  (2000 -2200  A)  where  the  Herz- 
berg  continuum  absorption  has  become  quite  weak,  and  the  ozone  absorption 
cross  section  at  the  short  wave  portion  of  the  Hartley  band  has  its  mini¬ 
mum  value.  In  this  spectral  region,  only  5  percent  of  the  solar  energy 
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Fig.  6.  The  spectral  distribution  of  depth  of  penetration  for  solar 

radiation  assuming  atmospheric  transmissivity  of  5  percent. 
The  two  curves  are  for  mean  conditions  corresponding  to 
30"  summer  and  60"  winter. 
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received  at  the  top  of  the  atmosphere  penetrates  to  levels  much  below 
35  km.  .At  the  time  of  high  sun,  the  solar  radiation  is  able  to  penetrate 
somewhat  deeper  into  the  stratosphei’e.  Although  3 -body  recombination 
at  these  heights  is  rather  rapid,  there  is  insufficient  energy  for  large- 
scale  dissociation  of  molecular  oxygen.  The  result  is  that  the  maximum 
equilibrium  amount  of  ozone  is  found  at  levels  just  under  35  km.  Below 
the  level  of  maximum  ozone  there  is  little  production  of  atomic  oxygen 
and  the  ozone  concentration  falls  rapidly. 

The  influence  of  the  strong  Hartley  absorption  bands  can  be 

O 

seen  in  the  shape  of  the  penetration  curve  in  the  region  2200-3000  A.  Very 

O 

little  solar  energy  of  wavelength  2550  A  can  penetrate  below  about  40  km. 
The  strong  absorption  even  at  the  top  of  the  stratosphere  (about  50  km) 
results  in  the  relatively  high  temperatures  found  at  the  stratopause. 

Again,  the  relative  height  difference  of  constant  atmospheric  trans  - 
missivity  in  this  spectral  interval  is  approximately  what  one  finds  for 
the  height  difference  of  the  level  of  maximum  temperature  between  equa¬ 
torial  regions  and  high  latitude  winter.  Since  we  have  seen  that  the  ozone 
concentration  is  not  sensitive  to  the  temperature  at  these  levels,  the 
temperature  structure  in  the  region  of  the  stratopause  must  result  from 
the  ozone  distribution  and  the  absorption  of  solar  energy  at  these  heights 
and  latitudes. 

The  two  curves  in  Fig.  6  cross  at  about  3000  A.  This  results 
from  the  fact  that  the  photochemical  distribution  rather  than  the  actual 
ozone  distribution  was  used  for  levels  below  25  km.  The  large  slope 
of  the  two  curves  near  3000 A,  however,  indicate  the  important  ab¬ 
sorption  transition  region  from  the  Hartley  to  the  Huggins  bands. 
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5.  Heating  of  the  upper  atmosphere 

The  processes  involved  in  the  reactions  described  in  Eqs  (1  -5) 
are  such  that  the  energy  exchanged  in  photochemical  dissociation  and  re¬ 
combinations  are  always  almost  in  balance  so  that  the  solar  energy  ab¬ 
sorbed  during  the  day  is  almost  completely  used  to  increase  the  kinetic 
energy  of  the  absorbing  molectiles.  This  results  in  a  temperature 
change  of  the  air  at  the  levels  v/here  the  energy  is  absorbed. 

The  total  absorption  per  unit  volume  at  any  level  by  molecular 
oxygen  and  ozone  is  then  given  by 


'2X^2  ■  °‘3X^3l‘^^ 

(15a) 

'2x^2  ■ 

(15b) 

where  represents  the  available  solar  energy  at  the  top  of  the 

d 

atmosphere  in  ergs  cm"®  sec'^A"’’ ,  and  the  integration  is  taken  over  all 
wavelengths  for  which  the  factors  signi- 

9  9 

ficantly  different  from  zero  (i.  e.  ,  1700  A  ^  X  ^2424  A  for  E2.  and 
2000  A  ^X^STBOGA  for  E^). 

The  heating  of  the  atmosphere  at  any  level  is  given  by  the  ab¬ 
sorption  per  unit  mass.  Thus,  the  total  heating  during  the  day  is 

(I) 


db> 


(16) 


where  u  is  the  hour  angle  of  the  sun  represents  the  local  hour 
angle  at  sunrise  or  sunset)  and  Ap  is  the  pressure  difference  between 
the  top  and  bottom  of  the  layer  in  which  the  absorption  occurs.  For  the 
computations  it  was  assumed  that  the  molecular  oxygen  concentration 
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is  constant  during  the  day,  that  the  local  ozone  concentration  is  also  con¬ 
stant  during  the  day  and  is  given  by  the  photochemical  equilibrium  dis¬ 
tribution  at  all  levels.  This  latter  assumption  is  probably  reasonably 
correct  down  to  about  40  km  below  v/hich  the  distribution  begins  to  de¬ 
part  radically  from  photochemical  equilibrium. 

Heating  computations  were  made  separately  for  molecular  exygen 
and  ozone.  It  was  found  that  atomspheric  heating  due  to  solar  energy 
absorption  by  molecular  oxygen  is  completely  negligible  in  the  strato¬ 
sphere  and  mesosphere  when  compared  to  heating  due  to  absorption  by 
ozone.  For  example,  the  heating  rate  due  to  oxygen  absorption  is  at 
most  1  -2  degrees  per  day  during  the  summer  at  about  80  km  and  de¬ 
creases  rapidly  with  decreasing  height.  On  the  other  hand,  heating  due 
to  ozone  absorption  is  of  the  order  of  10-20  degrees  per  day  through 
most  of  the  stratosphere  and  mesosphere  (except  for  high  latitudes  dur¬ 
ing  the  winter).  At  levels  above  80  km,  however,  the  oxygen  absorption 
becomes  important  and  as  the  ozone  concentration  goes  to  zero,  the 
heating  of  the  lower  layers  of  the  thermosphere  is  due  primarily  to  ab- 
sorption  by  molecular  oxygen  in  the  spectra  intervals  \  <  1800  A. 

The  distribution  with  height  and  latitude  of  the  calculated  heat¬ 
ing  rates  are  shown  as  a  summer-winter  cross  section  in  Fig.  7.  The 
values  are  given  in  *C/day  for  combined  absorption  due  to  molecular 
oxygen  and  ozone. 

It  can  be  seen  from  Fig.  7  that  the  rate  of  heating  is  a  maximum 
at  the  stratopause  and  is  about  Ib'C/day  at  high  latitudes  during  the 
summer  decreasing  to  about  3"C/day  at  60*  during  the  winter.  The 
level  of  maximum  heating  is  lowest  at  the  equator  (about  45  km)  and 
highest  (about  55  km)  at  high  latitudes  during  the  winter.  Above  this 
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Fig.  7.  The  distribution  of  heating  due  to  absorption  of  solar 
radiation  by  O3  and  Og  (“C/day). 
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region  of  maximum  heating  there  seems  to  be  a  shallow  layer  of  de¬ 
creased  heating  rate  and  then  a  deep  layer  of  strong  heating,  particularly 
at  high  latitudes  during  the  summer.  Below  45  km  the  heating  values 
are  everywhere  decreased  being  less  than  2°C/day  below  Z7  km. 

It  should  be  pointed  out  that  the  computed  heating  rates  are  some¬ 
what  unreliable  at  levels  below  about  35  km  because  of  the  str.jmg  de¬ 
partures  from  photochemical  equilibrium  concentrations. 

At  levels  above  60  km,  the  indicated  temperature  changes  due  to 
absorption  of  solar  radiation  could  be  incorrect  by  as  much  as  an  order 
of  magnitude.  At  these  elevations  the  ozone  distribution  again  departs 
somewhat  from  photochemical  equilibrium,  and  as  was  pointed  out  above, 
the  value  used  for  the  recombination  rate  ratio,  k^^/kg,  was  probably 
much  too  small.  In  addition,  secondary  reactions  involving  H,  OH,  and 
NO  have  been  neglected  in  the  photochemical  computations  (see,  for  in¬ 
stance,  Horiuchi,  1961;  Wallace,  1962).  These  effects  would  all  tend  to 
reduce  the  amount  of  ozone  present  and  reduce  the  heating  in  the  meso¬ 
sphere,  particularly  above  65  km.  The  computed  heating  rates  in  the 
region  of  35-60  km,  however,  are  quite  realistic  and  are  quite  clearly 
reflected  in  the  temperature  pattern  at  these  levels. 
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PART  III 


The  Distribution  of  Total  Ozone  in  the  Northern  Hemisphere 


by 

J.  London 


1 .  Introduction 

Measurements  of  total  ozone  in  a  vertical  column  in  the  atmo¬ 
sphere  have  been  made  since  about  1921  (Fowle,  1928).  The  first  set 
of  consistent  observations,  however,  resulted  from  a  standardized 
program  developed  in  1924  by  Professor  G.  M.  B.  Dobson  at  Oxford, 
England  (Dobson  and  Harrison,  1926).  This  observational  program 
was  based  on  a  spectrographic  technique  which  made  use  of  a  method 
first  suggested  by  Fabry  and  Buisson  (1921)  whereby  the  total  ozone 
was  calculated  from  observations  of  solar  radiation  (or  radiation  from 
any  stellar  source)  in  the  region  of  the  Huggins  ozone  absorption  bands 
(fti  3200  A).  Dobson  made  use  of  the  fact  that  if  the  received  light  in¬ 
tensity  is  measured  at  two  different  wavelengths  (one  of  strong,  the 
other  of  weak  absorption),  the  calculated  total  ozone  is  proportional 
to  the  logarithm  of  the  intensity  ratio  and  inversely  proportional  to 
the  difference  in  absorption  coefficients  at  the  two  different  wavelengths. 

In  1926  the  ozone  observing  program  was  extended  to  a  few  select¬ 
ed  places  in  both  hemispheres,  and  for  the  next  twenty -five  years  routine 
observations  were  taken  at  only  a  few  stations.  An  organized  obser¬ 
vational  program  involving  a  network  of  ozone  observing  stations 
(mostly  in  Europe)  was  started  by  the  then  International  Meteorological 
Organization,  in  1951.  At  the  start  of  the  IGY  program  (1957)  there 
were  about  twenty  stations  throughout  the  worid  taking  regular  daijy 
total  ozone  ohserva.tiona .  Since  then  this  niimber  has  more  than 
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doubled.  There  are  now  approximately  50  stations  observing  and  re¬ 
porting  total  ozone  values  on  a  routine  basis,  and  the  number  is  con¬ 
stantly  increasing. 

The  fundamental  observing  method  developed  by  Dobson  is  still 
used  today  by  the  extensive  ozone  observing  network,  although  a  number 
of  basic  improvements  in  the  instrumentation  were  made  with  the  in¬ 
troduction  in  1931  of  the  spectrophotometer  and  the  new  photo-electric 
instrument  in  1939.  The  presently  used  instrument  and  a  description 
of  the  technique  of  its  use  is  given  by  Dobson  (1957). 

Soon  after  the  first  observations  it  was  discovered  (Dobson, 
al.  1926;  Fowle,  1928,  1929)  that  the  vertically  integrated  total  ozone 
amount  had  large  day-to-day  variations  but  on  the  average  was  a  maxi¬ 
mum  during  the  spring  and  a  minimum  during  the  fall.  As  a  resxilt  of 
the  carefully  planned  observational  program,  Dobson  and  his  colleagues 
(Dobson,  et  al.  1928)  were  also  able  to  show  that  this  seasonal  variation 
prevailed  in  the  Southern  Hemisphere  and  that  total  ozone  values  were, 
on  the  average,  greatest  at  sub-polar  latitudes  (s»60")  and  least  near 
the  equator.  Their  studies  clearly  indicated  that  total  ozone  is  related 
to  synoptic  pressure  patterns  at  the  surface  and  aloft.  High  ozone  amounts 
tend  to  occur  with  low  surface  pressure;  low  ozone  amounts  tend  to  occur 
with  high  surface  pressure.  Later  studies  have  substantiated  these  re¬ 
lationships  for  various  regions  and  conditions.  Moreover,  it  has  been 
shown  that  the  correlations  between  ozone  and  pressure  and  temperature 
patterns  aloft  are  invariably  larger  during  winter  and  spring  than  during 
summer  and  fall  seasons. 

Although  the  strong  latitudinal  gradient  of  ozone  was  noticed 
early  in  ozone  research,  it  was  first  .suggested  by  Haurwitz  (1938) 


in  comparing  ozone  observations  taken  at  Shanghai  with  those  taken 
in  Western  Europe  that  there  might  be  a  longitude  variation  in  total 
ozone  as  related  to  surface  pressure  patterns.  This  possible  longi¬ 
tude  variation  was  again  noted  by  Craig  (1950)  and  has  since  been 
verified  by  many  for  specific  regions,  (See,  for  instance,  Kulkarni, 
et  al.  1959.) 

2.  Ozone  data 

The  increased  number  of  ozone  observations  in  recent  years 
has  made  it  possible  to  study  the  distribution  of  total  ozone  over  the 
entire  Northern  Hemisphere.  In  order  to  get  as  widespread  coverage 
of  the  data  as  possible  it  is  necessary  to  use  all  of  the  available 
ozone  observations  that  could  be  put  into  a  consistent  set. 

It  was  mentioned  earlier  that  the  calculated  ozone  amount 
is  inversely  proportional  to  the  adopted  value  of  the  difference  in 
ozone  absorption  coefficients  for  the  two  wavelengths  at  which 
measurements  are  taken.  That  is,  X  cm  of  ozone  is  proportional 
to  1  /(a  -  a'). 

The  early  observations  directed  by  Dobson  made  use  of 
absorption  coefficients  reported  by  Fabry  and  Buisson  (Dobson  et  al, 
1926).  Later,  results  of  Ny  and  Choong  (1932)  were  used  which  gave 
slightly  lower  computed  ozone  amounts.  More  recently,,  careful 
laboratory  measurements  by  Vigroux  (1953)  and  Inn  and  Tanaka  (1953) 
have  given  absorption  values  substantially  lower  than  those  reported 
by  Ny  and  Choong,  The  absorption  coefficients  of  Vigroux  have  been 
adopted  by  the  World  Meteorological  Organization  for  all  measurements 
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made  after  1  July  1957,  Consequently  all  observations  prior  to  that 
time  have  been  corrected  to  be  consistent  with  the  Vigroux  coefficients. 

Total  ozone  values  are  reported  here  in  Dobson  units 

-3 

(10  cm  NTP).  This  is  the  equivalent  depth  of  a  layer  of  ozone  if 
all  the  ozone  molecules  in  a  unit  vertical  column  were  brought  to 
normal  temperature  and  pressure  (0°C,  1013  mb). 

Although  there  are  a  few  partial  sources  for  total  ozone 
data,  (see  for  instance, Craig,  1950)  most  of  the  data  are  scattered 
throughout  the  literature,  much  of  it  still  unpublished.  An  attempt  was 
therefore  made  to  collect  and  tabulate  all  of  the  available  data  together 
with  the  sources  of  information.  In  this  endeavor  we  were  helped 
greatly  by  the  kind  cooperation  of  Sir  Charles  Normand  and  Professor 
G.  M.  B.  Dobson  then  secretary  and  president  repsecti  vely  of  the  In¬ 
ternational  Ozone  Commission.  The  sources  and  data  are  tabulated 
as  appendix  A  and  B  at  the  end  of  this  report.  The  source  list  is  given 
only  for  data  prior  to  the  start  of  the  ICY  ( 1  July  1957)  since  almost 
all  data  taken  subsequent  to  this  date  (July  1957  -  December  1959) 
have  been  tabulated  by  the  World  Meteorological  Organization.  Ar¬ 
rangements  are  currently  being  made  by  WMO  to  have  the  data  pub¬ 
lished  as  annual  summaries.  All  pre-IGY  data,  which  were  formally 
tabulated  according  to  the  Ny-Choong  absorption  values,  have  been 
multiplied  by  the  factor  1.  35  to  make  the  data  consistent  with  the  new 
Vigroux  scale  (see  Perl  and  Diiiisch,  1959). 

5.  The  hemispheric  ozone  distribution 

The  data  for  the  distribution  of  total  ozone  over  the  Northern 
Hemisphere  have  been  analyzed  for  the  whole  year  and  for  each  of  the 
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four  seasons.  The  data  used  were  the  totality  of  observations  made  at 
the  various  stations  from  the  start  of  their  observations  through  19  59, 
corrected  whan  necessary  to  be  consistent,  with  the  new  (Vigroux) 
absorption  coefficients  presently  being  used  with  the  Dobson  spectro¬ 
photometer  as  discussed  above.  The  analysis  is  given  in  Figs.  1(a)  - 
(d)  for  winter  (Dec,  Jan,  Feb);  spring  (Mar,  Apr,  May);  summer  (Jun, 
Jul,  Aug. ):  fall  (Sep,  Oct,  Nov);  and  in  Fig.  2  for  a  composite  annual 
distribution. 

Total  ozone  is  seen  to  be  a  minimum  at  the  equator  (about 
240  DU)  and  generally  to  increase  poleward  during  all  seasons.  In 
general,  the  latitudinal  gradient  is  a  maximum  during  the  spring  and 
minimum  during  the  fall.  In  polar  regions  there  are  three  centers  of 
high  ozone  amounts  found  at  about  60*N  over  Eastern  Canada,  North 
Central  Europe  and  Eastern  Siberia.  The  strength  of  these  centers 
increases  from  winter  to  spring  (to  a  maximum  of  about  460  DU)  and 
then  decreases  from  spring  to  fall.  Near  the  equator,  and  during  all 
seasons,  there  is  relatively  little  longitude  variation  in  total  ozone. 

For  instance,  Marcus  Island  (24*N  154°E)  and  Mavina  Loa  ( 19‘’N  155*W) 
have  practically  identical  ozone  amounts  during  all  seasons  of  the  year. 
At  higher  latitudes,  however,  this  is  not  so.  There  seems  to  be  at 
least  three  well  developed  ozone  ridges  extending  from  the  high  ozone 
centers  mentioned  above.  These  ridges  are  found  in  the  vicinity  of  the 
east  coa.sts  of  North  America  and  Asia  and  in  Central  Europe.  Although 
the  ridges  are  semipermanent  features  of  the  ozone  distribution,  they 
are  best  developed  during  winter  and  spring  [Figs.  1(a),  (b)]  and  least 
well  developed  during  summer  and  fall  [Figs.  1(c),  (d)  ]. 
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Comparison  between  the  maps  showing  hemispheric  distri¬ 
bution  of  ozone  and  maps  of  the  general  circulation  of  the  upper  tropo¬ 
sphere  and  lower  stratosphere  (Wahl  1958,  Panofsky  I960)  shows  that 
particularly  during  the  winter  and  spring,  there  are  three  well  developed 
contour  troughs  at  all  levels  from  200  mb  to  50  mb  located  in  the  same 
general  regions  as  for  the  ozone  ridges.  During  the  summer  and  fall, 
however,  although  the  contour  troughs  appear  at  200  mb,  they  have 
completely  disappeared  in  the  upper  level  maps.  There  is,  thus,  indica¬ 
tion  for  a  reasonably  strong  coupling  mechanism  between  the  development 
of  the  semipermanent  trough  pattern  at  upper  levels  and  the  total  ozone 
pattern  discussed  above.  The  dashed  lines  shown  in  Fig.  2  are  extra¬ 
polations  to  regions  where  there  are  no  data.  It  is  hoped  that  as  more 
data  become  available  the  area  of  uncertain  analysis  will  become  less 
and  less. 

The  low  ozone  anomaly  that  appears  over  the  Scandinavian 
peninsula  does  not  seem  to  have  any  well  defined  relation  to  the  upper 
air  circulation  pattern.  It  may  be  that  there  is  a  strong  local  bias  in 
the  observations  such  that  intervals  of  low  pressure  are  accompanied 
by  persistent  low  cloudiness  and  rain  which  would  prevent  accurate 
ozone  observations.  Because  of  the  inverse  relationship  between  pres¬ 
sure  and  ozone,  this  sampling  bias  would  result  in  unrepresentative  low 
ozone  values.  This  is  particularly  trvie  for  the  earlier  observations 
w/hen  the  various  cloud  corrections  were  not  clearly  defined.  More 
recent  data  (as  discussed  below  for  the  period  1958-59)  show  a  slight 
ozone  trough  in  the  general  vicinity  of  the  Scandinavian  peninsula  but 
not  nearly  as  pronounced  as  those  shown  in  Figs.  1(a)  -( J). 
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As  was  pointed  out  above,  in  an  attempt  to  get  widespread  cover¬ 
age  for  analysis  of  the  ozone  maps,  all  available  data  up  to  1959  were 
used.  In  this  grouping  of  data  there  were  many  stations  that  had  only 
one  or  two  years  of  observations  and  only  one  station  (Arosa,  Switzer¬ 
land)  that  had  over  30  years  of  data.  This  inhomogeneous  set  raises 
the  problem  of  the  representativeness  of  the  data.  The  two  longest  re¬ 
cords  of  data  are  those  for  Arosa  (32  years)  and  Tromso  (24  years). 
When  the  monthly  values  are  grouped  according  to  season,  the  number 
of  years  for  which  there  are  data  is  decreased.  The  least  number  of 
years  for  which  data  are  available  is,  of  course,  for  the  annual 
means.  The  following  table  gives  a  measure  of  the  inter-annual 
seasonal  and  annual  variation  of  ozone  amounts  for  three  stations  of 
longest  record. 


Table  1.  Standard  deviation  of  seasonal  and  annual 
total  ozone  values  (DU). 


Station 

Winter 

Spring 

Summer 

Fall 

Annual 

£ _ 

N 

(T  N 

cr 

N 

cr 

Ji 

cr 

N 

Tromso,  Norway(70“N) 

51 

12 

20  21 

15 

20 

15 

16 

18 

10 

Oxfort  ,  England  (52 'N) 

12 

9 

17  13 

10 

12 

12 

10 

8 

9 

Arosa,  Switzerland 
/47 

14 

29 

14  30 

8 

29 

7 

30 

9 

27 

It  should  be  pointed  out  that  in  most  cases  the  number  of  years 
available  to  compute  the  standard  deviation  was  quite  small  and  the 
values  given  in  Table  1  are  reliable,  if  at  all,  only  in  a  qualitative 
sense.  Three  general  characteristics  of  the  ozone  data,  however,  are 
suggested  by  the  numbers  in  Table  I.  The  variation  of  the  data; 

a)  increases  with  increasing  latitude 

b)  is  larger  during  winter  and  spring  than  summer  and  fall 

c)  is  about  20  DU  . 
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The  first  two  characteristics  mentioned  above  are  in  general 
agreement  with  our  knowledge  of  total  ozone  variations.  The  size  of  the 
standard  deviation  indicates  the  representativeness  of  the  analysis  shown 
in  Figs.  1(a) -(d)  and  Fig.  2.  Since  the  presence  of  the  high  ozone  cen¬ 
ters  and  the  approximate  placement  of  the  ozone  ridges  result  from  ozone 
differences  much  larger  than  20  DU  there  is  no  question  about  the  reality 
of  the  general  topography  of  the  longitude  variations.  However,  minor 
variations  in  the  ozone  distribution  probably  have  little  basis  in  fact. 

In  order  to  test  further  the  representativeness  of  the  distributions 
discussed  above,  a  separate  set  of  maps  was  analyzed  for  the  period 
1958-59.  These  are  shown  in  Figs.  3(a)  and  3(b)  for  the  spring  and 
fall.  Again  dashed  lines  indicate  areas  for  which  there  were  no  data. 
Although  the  data  available  for  analysis  of  Figs.  3(a),  (b)  were  consider¬ 
ably  less  than  for  Figs.  1, 2  the  advantages  are  that  the  data  were  co¬ 
herent  in  time  and  were  taken  with  improved  instruments  and  presumably 
more  accurate  techniques. 

The  results  shown  in  Figs.  3(a),  (b)  are  substantially  the  same 
as  discussed  above.  Again  there  is  a  pronounced  three  wave  pattern 
in  the  ozone  distribution  in  approximately  the  same  positions  as  shown 
before.  Also,  the  spring  patterns  are  much  stronger  than  those  for  fall. 
There  are  some  slight  differences  in  central  ozone  values  and  the  exact 
positions  of  the  ozone  ridges.  Although  it  is  probable  that  there  are  sig¬ 
nificant  year  to  year  variations  in  the  topography  of  the  ozone  distribution, 
the  present  set  of  observations  do  nox  permit  of  such  a  fine  analysis  as 
would  determine  such  variations.  The  differences  between  Figs.  1(b),  (d) 
and  Figs.  3(a),  (b)  are  certainly  much  smaller  than  their  similarities. 


Tlie  distribution  of  total  ozone  in  the  Northern  Hemisphere 

(1958-59). 


The  average  total  ozone  amount  around  each  latitude  circle  was 
computed  from  Figs.  l{a)-(d)  and  Fig.  2.  The  results  are  shown  in 
Fig.  4,  giving  mean  total  ozone  as  a  function  of  latitude  for  each  season 
and  for  the  annual  distribution.  The  data  did  not  extend  much  beyond 
20°  latitude  and  the  lower  parts  of  the  cvrrves  represent  extrapolated 
values.  Also,  the  fact  that  the  annual  curve  does  not  represent  the  aver¬ 
age  of  the  other  four  curves  in  the  region  north  of  75°  indicates  the  slight 
lack  of  consistency  of  the  analysis  in  this  region.  In  general,  however, 
the  patterns  hold  together  quite  well. 

During  all  seasons  there  is  increasing  ozone  to  the  north  with  the 
maximum  indicated  at  about  60 “-70°.  The  largest  gradient  is  found  in 
middle  latitudes  (30  “-50“)  during  all  seasons  but  predominantly  during 
spring.  Also,  it  is  most  noticeable  that  the  seasonal  variation  is  small¬ 
est  in  equatorial  latitudes  and  largest  in  polar  latitudes.  If  the  ozone 
amounts  are  averaged  over  the  hemisphere  for  each  of  the  seasons  we 
find  that  the  observed  distribution  gives: 


Winter 

298  DU 

Spring 

330  DU 

Summer 

297  DU 

Fall 

274  DU 

Annual 

299  DU 

This  can  be  compared  to  the  amount  of  ozone  which  would  be  maintained 

in  a  photochemical  equilibrium  atmosphere.  If  we  integrate  the  vertical 

distribution  of  equilibrium  ozone  as  given  in  the  preceding  section  for 

each  latitude  and  the  two  seasons  we  get 

Summer  -  347  DU 

Winter  -  171  DU 

Annual  -  259  DU 


TOTAL  OZONE  -  D.U.  ( 10  cm.  STP) 


80 


440 


400 


360 


320 


280 


240 
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Although  the  annual  amount  is  only  slightly  less  than  the  observed  annual 
amount,  of  most  concern  is  the  marked  differences  in  seasonal  and  lati¬ 
tudinal  variations.  Thus,  the  observable  amounts  of  ozone  can  be  pro¬ 
duced  by  photochemical  processes  but  the  temporal  and  latitudinal  dis  - 
tributions  most  certainly  are  controlled  by  the  large  scale  general  circu¬ 
lation  of  the  upper  troposphere  and  lower  stratosphere.  The  large  scale 
distribution  of  total  ozone  shows  a  marked  relation  to  the  large  scale 
atmospheric  pressure  systems.  Also  ozone  must  be  transported  from 
the  region  of  production  (equatorial  latitudes)  to  polar  regions  to  account 
for  the  observed  ozone  distribution.  Both  of  these  facts  indicate  that  it 
is  most  likely  that  the  atmospheric  transporting  mechanism  for  ozone  is 
essentially  the  same  as  that  responsible  for  maintaining  the  heat  and 
momentum  budgets  of  the  atmosphere.  It  is  of  considerable  general 
interest  therefore,  to  probe  deeper,  as  more  data  become  available, 
into  the  problem  of  the  relation  of  the  general  circulation  of  the  atmo¬ 
sphere  and  the  variations  in  the  ozone  budget. 
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Appendix  A 

Preliminary  Data  Source  -  Total  Ozone 
(Northern  Hemisphere) 

Data  Prior  to  1  July  1957 


Station 

Lat. 

Long. 

Period  Reference 

Aarhus,  Denmark 

5618N 

1037E 

Sep  51 -Nov  53 

18 

Mar  54-Jun  57 

18 

Abisko,  Sweden 

6820N 

1847E 

Jul  26,  Apr-Sep  27 

10 

Dec  34,  Jan-Mar  35 

2 

Ahmedabad,  India 

2304N 

7238E 

Jan-Apr  55 

13 

Albuquerque, 

New  Mexico 

3503N 

10637W 

Jan-Jul  54 

22 

Aldergrove 

North  Ireland 

5439N 

0613W 

Jan  52 -Oct  56 

18 

Jan  -  Apr  57 

18 

Arosa,  Switzerland 

4647N 

0941E 

Jul  26-Jun  29 

19 

Jul -Aug  30 

19 

Aug  31 -Sep,  Nov  53 

19 

Mar  54-Jun  57 

19 

Bismarck, 

North  Dakota 

4646N 

10045W 

Jan  54-Nov  56 

22 

Bombay,  India 

1856N 

7254E 

Oct  36-Sep  38 

6 

Cagliari,  Sardinia 

3914N 

0906E 

Oct  54-Nov  55 

18 

Jan-Feb  56 

18 

Jul  56-Jun  57 

18 

Camborne,  England 

5013N 

0519W 

Nov  50-May  56 

18 

Oct  56«Jun  57 

18 

College,  Alaska 

6444N 

14730W  Apr-Oct  52 

5 

Feb-Oct  53 

5 

86 


Station 

Lat. 

Long. 

Period  Referenc 

Dombas,  Norway 

6208N 

0908E 

Mar -Nov  40 

16 

Feb-Nov  41 

16 

Mar-Sep  42 

16 

Jan-Dec  43 

16 

Mar  44-Jun  46 

16 

Downham  Market 

England 

5236N 

0023E 

Nov  50-Nov  5 1 

18 

Edmonton,  Canada 

5334N 

11331W 

Aug  50-Dec  52 

12 

Eskdalemiur 

Scotland 

5519N 

0312W 

Jun  57 

18 

Flagstaff,  Arizona 

3508N 

11140W 

Sep  54-Dec  56 

1 

Fort  Worth,  Texas 

3250N 

9703W 

Jan  54-Nov  55 

22 

Gulmarg,  India 

3403N 

7224E 

Jun-Dec  55 

13 

Apr-Sep  56 

13 

Helwan,  Egypt 

2952N 

3118E 

Sep  2 8 -May  29 

9 

Jul-Oct  29 

9 

Hemsby,  England 

5241N 

0141E 

Nov-51 -Aug  52 

18 

Nov  52»Sep  55 

18 

Kabul,  Afghanistan 

3432N 

6910E 

Jan-Dec  51 

15 

Kodaikanal,  Indi.a 

1014N 

7728E 

Sep  28-Aug  29 

9 

Lerwick,  Scotland 

6008N 

OlllW 

Jul-Sep  26 

7 

Mar-Oct  27 

7 

May -Aug  51 

18 

Mar  52-Nov  54 

13 

Jan  55-Jun  57 

18 

Lindenberg,  Germany  52  12N 

1408E 

Jul-Oct  26 

10 

Apr -Oct  27 

10 

Magny  Les  Hameaux, 

France 

4844N 

0204E 

Jan  55-Jun  57 

18 

Marseille,  France 

4317N 

0522E 

Jan-Dec  27 

3 

Jan-Dec  28 

4 

Messina,  Sicily 

3812N 

1533E 

Jul  54-Jun  57 

18 

Mount  Abu,  India 

2436N 

7243E 

Jan  55-Jun  57 

13 

Station 

Lat. 

LiOng. 

Period  Reference 

New  DelJii. India 

2838N 

7713E 

Nov  45-Mar  47 

14 

Jan  55-Jun  87 

13 

New  York,  N.  Y. 

4046N 

7348W 

Jan  41 -Nov  42 

8 

Jan-Jul  43 

8 

Oct  43 -Dec  44 

8 

Oslo,  Norway 

6005N 

lOOOE 

Jul  46 -Apr  49 

16 

Oxford,  England 

5146N 

0116W 

Feb-Nov  25 

7 

F eb-Oct  26 

10 

Mar-Oct  27 

10 

Mar  -Oct  28 

10 

Nov  50-Jun  57 

18 

Pasadena,  Calif. 

3409N 

11810W 

Sep-Dec  39 
Feb-Mar,  May- 

7 

Aug  40 

7 

Potsdam,  Germany 

5223N 

1304E 

Feb-Apr,  Jun  42 

Jan-Mar,  May,  Aug 

Sep  43 

Reykjavik,  Iceland 

6408N 

2154W 

Feb-Oct  52 

18 

Feb-Oct  53 

18 

Feb-Oct  54 

18 

Mar -Jul  55 

18 

Rome,  Italy 

4154N 

1228E 

Apr  54-Jun  57 

18 

Santa  Maria,  Arores 

3655N 

2510W 

Jun,  Jul  51 

18 

Jan  53-Jul  56 

18 

Shanghai,  China 

3111N 

12129E 

Jan-32-Dec  33 

17a 

Jan  34-Dec  36 

17b 

Jan  -Dec  37 

17c 

Feb  -SS-Dec  42 

17c 

Spitzbergen,  Norway  7800N 

1500E 

Jul-Aug  29 

7 

Nov  50-Sep  51 

18 

Dec  51 -Aug  53 

18 

Nov  63-Oct  54 

18 

Dec  54 

18 

Mar -Sep  55 

18 

Mar -Aug  56 

18 

Srinagar,  India 

3405N 

7450E 

Jan-Apr  56 

13 

Nov  56~Jun  57 

13 

88 


Stahi  on 

Lat. 

Long. 

Period 

Reference 

Table  Mountain, 

California 

3424N 

11742W 

Aug  2 8 -Aug  29 

9 

Tateno,  Japan 

3603N 

14008E 

Jul  55-Nov  56 

20 

Feb-Jun  57 

20 

Tromso,  Norway 

6940N 

1857E 

Jul  35  -Oct  36 

21 

Jan-Dec  37 

21 

Mar-Oct  38 

21 

Mar  39 "Dec  42 

21 

Jan  43 -Nov  46 

16 

Jan  47-Nov  48 

16 

Jan,  Mar-May  49 

16 

Mar  51,  Dec  52 

18 

Apr-Oct  53 

18 

Feb-Oct  54 

18 

Feb-Oct  55 

18 

Feb-Jun  56 

18 

Uccle,  Belgium 

5048N 

0421E 

Sep  51 -Aug  52 

18 

Mar  54-Oct  55 

18 

Apr-Jun  57 

18 

Uppsala,  Sweden 

5952N 

1737E 

Jul-Oct  51 

18 

Mar  52-Nov  54 

18 

Jan  55-Jun  57 

18 

Valentia, Ireland 

5156N 

1015W 

Sep-Oct  26 

10 

Mar-Oct  27 

10 

Washington,  D.  C. 

3851N 

7702W 

Feb  48-Mar  50 

11 

Jan  54»Sep  56 

22 

Data  Source  Reference  List 


Total  Ozone 
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APPENDIX  B 

Monthly  means  of  observed  total  ozone  (Dobson  units  -  10  ^  cm  NTP). 
All  values  have  been  converted  to  be  consistent  with  Vigroux  absorption 
coefficients.  Number  of  days  of  observations  for  each  month  is 
indicated  when  available. 
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Aarhus,  Denmark,  56*18'N,  10°37'E 


Year 

Jan 

Feb 

Mai 

•  Apr 

May 

Jun 

Jul 

Ai^ 

Se£ 

Oct 

Nov 

Dec 

21 

30 

~Z4 

31 

1951 

287 

264 

269 

314 

28 

27 

23 

29 

28 

28 

27 

26 

26 

20 

29 

27 

1952 

369 

388 

443 

424 

416 

407 

352 

318 

321 

332 

309 

362 

28 

21 

31 

24 

25 

26 

20 

27 

18 

27 

21 

1953 

344 

386 

382 

390 

405 

389 

378 

339 

292 

296 

283 

17 

28 

30 

27 

28 

28 

26 

23 

22 

26 

1954 

416 

454 

419 

393 

392 

366 

340 

303 

298 

306 

29 

27 

29 

28 

28 

27 

30 

31 

23 

25 

22 

24 

1955 

352 

446 

431 

411 

420 

397 

356 

359 

313 

306 

279 

321 

26 

27 

27 

24 

19 

28 

30 

25 

28 

25 

25 

28 

1956 

350 

426 

397 

393 

420 

392 

360 

374 

313 

305 

303 

290 

24 

22 

27 

29 

28 

20 

1957 

326 

445 

397 

408 

426 

378 

Abisko,  Sweden,  68"20'N,  ISMT'E 

8 

1926  301 

2  14  7  9  5  7 

1927  449  424  347  321  313  286 

1928 

1929 

1930 

1931 

1932 

1933 

3 

1934  237 

11  12  6 

1935  296  363  356 


Ahmedabad,  India,  23  *04'N,72 ‘38'E 

28  26  30  29 

1955  222  227  235  242 


93 


Albuquerque,  New  Mexico,  35‘’03'N,  106*37'W 


Year 

Jan 

Feb 

Mar 

Miy 

Jun 

Jul 

Sep. 

Oct 

Nov 

Dec 

17 

22 

14 

4 

16 

17 

3 

1954 

307 

328 

352 

354 

362 

332 

326 

Aide 

rgrove 

,  N.  Ireland.  54 

“39'N, 

06* 

13'W 

4 

26 

29 

30 

18 

21 

31 

31 

29 

26 

28 

27 

1952 

400 

339 

416 

397 

385 

365 

335 

326 

303 

296 

267 

298 

28 

26 

31 

30 

31 

29 

29 

13 

26 

31 

29 

25 

1953 

286 

324 

343 

375 

364 

362 

359 

343 

288 

282 

260 

295 

27 

21 

31 

30 

30 

30 

30 

30 

30 

29 

27 

22 

1954 

330 

375 

367 

407 

407 

390 

370 

345 

313 

287 

280 

276 

23 

23 

30 

28 

31 

29 

31 

30 

30 

31 

28 

19 

1955 

325 

416 

377 

358 

397 

381 

351 

328 

291 

287 

261 

305 

18 

25 

21 

26 

28 

24 

24 

23 

17 

7 

1956 

335 

348 

384 

394 

396 

377 

369 

374 

314 

290 

12 

14 

19 

9 

1957 

314 

409 

397 

381 
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Arosa,  Switzerland,  46*47'N,  09°41'E 


Y  ear 

Jan 

Feb 

Mar 

May 

Jun 

Jul 

Se£ 

Oct 

Nov 

Dec 

5 

28 

27 

25 

20 

13 

1926 

303 

292 

275 

268 

286 

313 

19 

24 

24 

18 

31 

27 

30 

29 

26 

29 

22 

23 

1927 

351 

351 

362 

360 

367 

337 

326 

298 

294 

273 

275 

305 

19 

21 

17 

21 

19 

28 

29 

26 

19 

21 

17 

22 

1928 

313 

313 

335 

347 

360 

318 

296 

286 

277 

258 

271 

302 

20 

25 

22 

25 

18 

6 

1929 

336 

386 

343 

401 

354 

333 

“ 

-  no  ' 

observations 

.... 

.... 

1 

11 

1930 

---no  observations  - 

— 

(279) 

286 

— 

no  obs. 

— 

7 

20 

19 

14 

21 

1931 

- no  observations  - 

— 

— 

301 

324 

273 

277 

291 

22 

25 

19 

11 

19 

22 

16 

24 

15 

16 

16 

12 

1932 

307 

,335 

356 

378 

347 

335 

322 

290 

283 

286 

283 

272 

14 

15 

24 

13 

16 

21 

27 

24 

21 

7 

11 

21 

1933 

344 

350 

3  84 

366 

374 

359 

322 

307 

298 

302 

301 

325 

12 

20 

19 

23 

21 

17 

24 

19 

21 

17 

19 

12 

1934 

322 

310 

377 

345 

350 

343 

317 

310 

272 

277 

282 

287 

10 

12 

21 

18 

23 

20 

17 

8 

14 

11 

12 

14 

1935 

320 

374 

354 

369 

360 

309 

320 

301 

279 

267 

290 

302 

12 

18 

19 

14 

15 

18 

19 

24 

20 

21 

16 

23 

1936 

317 

371 

382 

369 

359 

339 

317 

305 

292 

299 

287 

296 

14 

10 

21 

15 

22 

19 

16 

15 

16 

8 

18 

15 

1937 

333 

339 

379 

367 

352 

337 

313 

311 

291 

273 

271 

343 

13 

21 

23 

21 

14 

15 

19 

8 

20 

20 

25 

17 

1938 

325 

355 

314 

377 

370 

324 

314 

314 

286 

270 

276 

288 

14 

19 

12 

11 

6 

16 

22 

12 

14 

17 

20 

17 

1939 

309 

329 

370 

335 

367 

326 

320 

302 

277 

294 

275 

299 

20 

13 

26 

24 

23 

21 

19 

21 

19 

10 

17 

21 

1940 

373 

385 

401 

412 

388 

373 

333 

311 

299 

282 

291 

325 

19 

21 

22 

21 

22 

26 

23 

27 

25 

21 

21 

20 

1941 

348 

379 

400 

393 

401 

347 

336 

324 

295 

290 

298 

298 

19 

20 

20 

18 

24 

21 

24 

24 

27 

19 

19 

25 

1942 

3  85 

405 

359 

392 

362 

335 

314 

298 

275 

264 

288 

302 

15 

15 

29 

26 

22 

22 

21 

2Z 

24 

27 

13 

18 

1943 

325 

329 

370 

350 

336 

339 

324 

292 

271 

283 

292 

309 

20 

16 

13 

21 

24 

25 

24 

26 

24 

24 

15 

25 

1944 

290 

351 

370 

347 

336 

337 

309 

295 

282 

284 

288 

309 

23 

20 

23 

26 

28 

29 

29 

28 

22 

27 

23 

9 

1945 

363 

347 

347 

359 

362 

339 

317 

316 

287 

279 

286 

302 

(Cotit.  ) 
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Arosa,  Switzerland  (Continued) 


Year 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

24 

18 

26 

25 

25 

22 

24 

16 

22 

28 

24 

16 

1946 

324 

339 

366 

348 

341 

330 

306 

296 

268 

282 

287 

306 

25 

19 

21 

25 

27 

24 

25 

30 

28 

28 

19 

14 

1947 

369 

382 

378 

355 

348 

329 

322 

313 

296 

287 

269 

301 

20 

20 

27 

26 

28 

23 

21 

27 

16 

24 

22 

25 

1948 

329 

356 

335 

350 

340 

333 

332 

301 

288 

272 

276 

310 

14 

26 

26 

27 

23 

29 

28 

23 

26 

27 

17 

20 

1949 

321 

351 

363 

344 

356 

341 

324 

310 

275 

264 

286 

283 

22 

20 

27 

23 

28 

27 

31 

24 

26 

25 

19 

17 

1950 

340 

351 

352 

367 

360 

341 

311 

305 

283 

279 

277 

328 

20 

20 

24 

29 

26 

26 

26 

26 

26 

30 

18 

24 

1951 

326 

386 

400 

382 

369 

351 

320 

310 

287 

288 

269 

301 

21 

22 

21 

23 

25 

27 

30 

27 

25 

19 

18 

17 

1952 

363 

369 

394 

371 

370 

345 

329 

309 

306 

287 

292 

321 

25 

25 

30 

26 

24 

22 

26 

22 

25 

3 

25 

no 

1953 

324 

362 

359 

367 

367 

345 

316 

306 

283 

(269) 

258 

obs 

23 

16 

24 

12 

25 

22 

12 

25 

24 

13 

1954 

-no 

oba  - 

-360 

400 

374 

348 

335 

317 

295 

276 

282 

275 

21 

11 

27 

24 

23 

26 

18 

28 

23 

24 

24 

11 

1955 

305 

360 

384 

359 

348 

339 

328 

321 

290 

«83 

269 

299 

17 

22 

21 

17 

21 

11 

23 

12 

23 

23 

21 

20 

1956 

325 

381 

363 

374 

348 

343 

315 

296 

277 

274 

275 

294 

21 

23 

29 

24 

19 

26 

1957 

325 

328 

338 

354 

367 

335 

Bismarck, 

,  N. 

Dak.  , 

46'46’NJOO 

”45' 

'W_ 

22 

22 

21 

19 

25 

21 

29 

30 

24 

26 

25 

25 

1954 

360 

363 

408 

416 

405 

363 

359 

359 

326 

307 

280 

321 

23 

18 

26 

18 

23 

25 

26 

28 

25 

29 

21 

23 

1955 

347 

373 

422 

368 

324 

370 

341 

316 

299 

272 

284 

288 

21 

19 

25 

25 

24 

26 

28 

28 

24 

30 

21 

1956 

295 

362 

375 

^94 

359 

345 

345 

329 

306 

282 

287 

96 


Bombay,  India.  18»55'N,  72°54'E 


Year 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Se£ 

Oct 

Nov 

Dec 

10 

2 

9 

1936 

260 

235 

256 

26 

17 

29 

28 

3 

2 

4 

14 

15 

28 

22 

31 

1937 

257 

280 

252 

262 

228 

226 

215 

233 

256 

262 

246 

248 

30 

28 

31 

27 

28 

9 

11 

10 

16 

1938 

246 

256 

257 

275 

268 

264 

267 

248 

256 

Cagliari,  Sardinia,  39 

“M'N. 

09*06iE 

31 

30 

31 

1954 

295 

299 

292 

30 

27 

28 

28 

31 

30 

31 

30 

30 

31 

27 

1955 

302 

321 

356 

370 

321 

324 

316 

314 

303 

292 

313 

9 

17 

31 

30 

30 

27 

27 

31 

1956 

313 

385 

314 

301 

299 

294 

301 

302 

29 

28 

31 

27 

30 

28 

1957 

310 

301 

321 

318 

366 

337 

Camborne, 

England,  ! 

50'13'N,  05 

*19'W 

16 

25 

1950 

292 

318 

29 

27 

26 

28 

31 

26 

28 

30 

28 

2 

5 

24 

1951 

329 

3  82 

385 

366 

409 

360 

320 

326 

282 

275 

271 

307 

21 

24 

26 

27 

28 

28 

30 

29 

28 

26 

23 

23 

1952 

354 

340 

397 

394 

423 

392 

360 

347 

335 

305 

286 

302 

20 

25 

30 

26 

26 

17 

27 

27 

25 

29 

27 

30 

1953 

301 

333 

366 

396 

384 

371 

355 

324 

283 

284 

264 

280 

27 

27 

28 

29 

29 

26 

29 

29 

27 

27 

24 

27 

1954 

320 

362 

366 

405 

394 

381 

355 

339 

301 

282 

277 

276 

28 

23 

29 

26 

24 

27 

29 

28 

26 

26 

29 

28 

1955 

317 

377 

378 

352 

379 

362 

344 

324 

282 

271 

257 

295 

25 

29 

29 

19 

8 

26 

25 

27 

1956 

326 

350 

350 

378 

348 

264 

262 

275 

27 

25 

28 

30 

30 

30 

1957 

301 

359 

325 

333 

350 

333 

97 


College,  Alaska  ,  64°44'N,  147°30»W 


Year 

Jan  Feb 

Mar 

Apr 

May 

Jun 

Se£ 

Oct  Nov  Dec 

9 

17 

19 

13 

16 

3 

6 

1952 

446 

411 

356 

326 

301 

299 

299 

5 

25 

27 

22 

18 

20 

6 

13 

11 

1953 

396 

438 

416 

394 

351 

287 

282 

282 

326 

10 

18 

15 

18 

25 

19 

22 

13 

8 

1954 

427 

437 

449 

411 

394 

370 

330 

296 

295 

Dombas,  Norv/ay, 

62* 

08'N, 

09“08'E 

9 

19 

22 

12 

10 

20 

24 

24 

6 

1940 

390 

392 

384 

322 

311 

305 

295 

264 

290 

6 

23 

28 

25 

23 

18 

25 

27 

22 

3 

1941 

416 

422 

422 

393 

324 

306 

318 

264 

287 

275 

15 

24 

23 

24 

30 

28 

10 

1942 

400 

364 

362 

347 

332 

301 

275 

13 

9 

19 

26 

30 

19 

27 

25 

26 

30 

26 

22 

1943 

290 

307 

373 

379 

369 

322 

314 

302 

273 

250 

269 

305 

4 

23 

27 

27 

22 

20 

30 

27 

23 

17 

1944 

363 

356 

358 

329 

302 

269 

283 

258 

269 

292 

7 

4 

24 

29 

22 

26 

29 

30 

29 

28 

24 

22 

1945 

311 

367 

367 

398 

386 

340 

305 

291 

277 

248 

253 

313 

29 

17 

27 

29 

31 

17 

1946 

318 

333 

386 

374 

254 

336 

Downham 

Market,  England,  52 

“36'N, 

00*23 

I'E 

26 

31 

1950 

302 

326 

28 

25 

31 

29 

31 

29 

31 

31 

29 

30 

17 

1951 

352 

420 

419 

420 

423 

394 

363 

350 

295 

291 

280 

Pdmonton,  Alberta, 

Can. 

,  5334'N, 

113° 

31  'W 

15 

16 

5 

9 

10 

1950 

310 

264 

283 

307 

295 

14 

16 

21 

18 

20 

18 

26 

23 

24 

23 

21 

12 

1951 

333 

343 

405 

359 

389 

401 

354 

330 

292 

290 

294 

322 

20 

25 

23 

21 

22 

18 

21 

20 

25 

29 

23 

21 

1952 

371 

398 

419 

398 

384 

350 

336 

321 

286 

268 

307 

299 

98 


Flagstaff.  Arizona.  35°08'N.  IllMO'W 


Year 

Jan 

Feb 

Mar 

Miy 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

22 

24 

17 

14 

1954 

268 

252 

257 

279 

11 

14 

19 

17 

16 

22 

23 

15 

24 

22 

13 

7 

1955 

275 

287 

290 

309 

306 

292 

279 

283 

267 

254 

252 

250 

7 

10 

19 

7 

20 

17 

10 

19 

16 

13 

6 

10 

1956 

282 

311 

306 

321 

306 

288 

279 

277 

267 

261 

245 

275 

Fort  Worth, 

Texas 

,  32 

"50'N, 

97003 

'W 

16 

26 

22 

19 

19 

27 

26 

27 

28 

17 

24 

21 

1954 

273 

294 

294 

303 

328 

316 

299 

283 

269 

253 

254 

265 

12 

18 

19 

20 

24 

26 

26 

25 

28 

27 

7 

1955 

276 

265 

269 

283 

298 

298 

292 

284 

269 

253 

267 

Gulmarg, 

India, 

34-03 

'N, 

72-24'E 

18 

21 

25 

28  26 

28 

24 

1955 

258 

243 

230 

233  235 

235 

249 

2  30 

24 

25 

22 

25 

1956 

280  257 

256 

250 

224 

222 

Helwan, 

Egypt,  29 

*52'N, 

31”: 

18'E 

3 

9 

16 

22 

1928 

256 

256 

248 

253 

23 

18 

15 

17 

2 

8 

20 

23 

6 

1929 

279 

283 

299 

294 

284 

275 

276 

277 

257 

Hemsby 

,  England, 

52° 

41 'N, 

01”41 

'E 

7 

23 

1951 

269 

314 

26 

27 

28 

29 

31 

30 

31 

12 

16 

30 

1952 

375 

360 

423 

415 

415 

418 

369 

354 

•' 

302 

314 

31 

27 

31 

30 

29 

30 

31 

30 

29 

24 

19 

27 

1953 

298 

348 

358 

382 

362 

375 

354 

316 

288 

280 

261 

296 

22 

20 

27 

30 

28 

30 

30 

27 

27 

22 

23 

8 

1954 

325 

397 

378 

412 

400 

373 

360 

339 

310 

280 

273 

299 

26 

20 

26 

25 

18 

24 

30 

30 

17 

1955 

321 

401 

374 

348 

384 

367 

344 

317 

284 

99 


Kabul,  Afghanistan,  34°32'N,  69“J0'E 


Year 

Jan 

Feb 

Mar 

May 

Jun 

Jul 

Sep 

Oct 

Nov 

Dec 

6 

10 

8 

10 

4 

5 

10 

18 

20 

13 

10 

11 

1951 

(273) 

268 

(314) 

321 

(308)  (236) 

249 

282 

260 

291 

271 

347 

Kodaikanal, 

,  India,  1 0 

'14'N, 

77” 

28'E 

17 

12 

7 

12 

1928 

24 

24 

23 

15 

18 

7 

14 

19 

246 

242 

231 

226 

1929 

230 

233 

238 

245 

253 

250 

253 

245 

Lerwick,  Scotland,  60 

"08'N, 

01  “ 

ll'W 

21 

17 

15 

1926 

328 

321 

301 

12 

13 

17 

14 

22 

20 

19 

5 

1927 

416 

416 

397 

341 

333 

335 

305 

307 

5 

26 

30 

18 

1951 

359 

362 

355 

336 

16 

30 

31 

30 

31 

31 

30 

31 

22 

3 

1952 

420 

404 

389 

398 

352 

337 

310 

302 

256 

250 

6 

27 

31 

30 

31 

30 

30 

31 

30 

30 

16 

4 

1953 

292 

337 

355 

393 

382 

366 

354 

326 

287 

273 

280 

246 

4 

27 

31 

29 

31 

30 

31 

29 

29 

30 

14 

1954 

318 

370 

385 

411 

385 

364 

352 

329 

310 

288 

284 

1 

27 

25 

30 

7 

2 

27 

31 

25 

30 

17 

2 

1955 

291 

404 

378 

355 

385 

371 

325 

306 

291 

273 

248 

260 

3 

28 

30 

30 

31 

30 

31 

31 

30 

31 

17 

1 

1956 

290 

313 

384 

409 

381 

378 

344 

348 

302 

277 

260 

248 

5 

23 

31 

30 

31 

30 

1957 

313 

438 

401 

381 

407 

369 

1926 


Lindenberg,  Germany,  5Z‘I2'N,  14°08'E 


23 

22 

19 

10 

320 

320 

294 

295 

7 

26 

22 

20 

19 

13 

7 

381 

403 

367 

335 

313 

316 

301 

1927 


100 


Magny  Les  Hameaux,  France,  48“44'Ni  02°04'E 


Year 

Jan 

Feb 

Mar 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

7 

15 

15 

20 

20 

19 

14 

7 

17 

17 

17 

11 

1955 

298 

354 

340 

335 

347 

336 

343 

297 

273 

271 

253 

273 

19 

18 

21 

12 

11 

19 

17 

12 

15 

17 

17 

11 

1956 

3  39 

369 

356 

381 

328 

337 

317 

340 

284 

288 

284 

295 

17 

14 

15 

16 

14 

5 

1957 

326 

351 

326 

350 

369 

366 

Marseille, 

France,  43 

•17'N, 

05“2 

2'E 

3 

17 

14 

29 

19 

24 

29 

27 

22 

24 

11 

7 

1927 

370 

389 

422 

432 

424 

378 

354 

330 

330 

321 

273 

316 

20 

22 

13 

17 

26 

24 

30 

14 

22 

21 

23 

19 

1928 

343 

378 

412 

405 

394 

351 

336 

313 

316 

313 

326 

354 

Messina, 

Italy, 

38* 

12'N, 

15‘33 

'E 

31 

31 

30 

31 

30 

31 

1954 

287 

284 

279 

271 

272 

275 

31 

28 

31 

30 

31 

30 

31 

31 

29 

30 

30 

31 

1955 

282 

280 

317 

335 

305 

298 

294 

284 

276 

267 

265 

276 

31 

27 

30 

28 

30 

25 

30 

31 

28 

29 

28 

30 

1956 

301 

354 

330 

341 

326 

340 

311 

292 

275 

265 

261 

303 

31 

28 

30 

30 

31 

30 

1957 

344 

332 

378 

397 

407 

351 

Mount  Abu, 

India 

,  24* 

36'N, 

72"43'E 

31 

27 

31 

30 

29 

27 

27 

11 

12 

29 

25 

30 

1955 

222 

227 

233 

241 

252 

246 

243 

243 

246 

242 

238 

231 

31 

29 

31 

27 

31 

24 

1 

9 

24 

11 

31 

1956 

237 

254 

254 

258 

256 

248 

246 

230 

230 

231 

224 

30 

27 

30 

6 

23 

1957 

219 

219 

242 

268 

248 

101 


New 

Delhi, 

India,  28 

“JS'N, 

77° 

13'E 

ISH. 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

19 

27 

1945 

230 

231 

28 

22 

28 

30 

29 

28 

21 

21 

28 

21 

29 

28 

1946 

231 

235 

249 

256 

254 

25? 

243 

235 

239 

243 

238 

248 

18 

26 

24 

1947 

261 

268 

269 

1948 

1949 

1950 

1951 

1952 

1953 

1954 

28 

27 

31 

30 

30 

22 

25 

25 

28 

26 

30 

31 

1955 

264 

257 

248 

257 

L'2 

253 

256 

250 

249 

257 

246 

256 

28 

29 

31 

30 

29 

26 

21 

25 

28 

27 

29 

31 

1.956 

275 

283 

277 

262 

264 

260 

269 

279 

265 

237 

234 

241 

25 

28 

30 

29 

31 

28 

1957 

250 

235 

242 

233 

243 

256 

New  York,  N,  Y.  , 

40“46'N,  73 

“46'W 

24 

27 

26 

28 

31 

26 

28 

31 

30 

30 

25 

30 

1941 

366 

417 

431 

388 

411 

390 

374 

358 

324 

308 

311 

325 

25 

22 

10 

17 

13 

4 

14 

24 

21 

5 

3 

1942 

3  84 

405 

420 

418 

388 

373 

354 

340 

307 

335 

325 

19 

8 

10 

16 

19 

21 

11 

357 

343 

11 

13 

13 

1943 

363 

386 

389 

400 

385 

333 

371 

msg 

msg 

330 

315 

339 

20 

23 

23 

17 

19 

26 

31 

24 

24 

26 

22 

23 

1944 

356 

384 

358 

408 

371 

362 

364 

348 

317 

311 

325 

324 

Oslo, 

Norway, 

60°05 

'  N,  10 

°00'E 

5 

3 

22 

26 

21 

23 

1946 

333 

295 

276 

241 

261 

290 

26 

24 

26 

27 

27 

25 

27 

26 

29 

26 

25 

26 

1947 

284 

431 

430 

388 

350 

343 

326 

282 

282 

272 

276 

337 

27 

25 

27 

30 

30 

28 

27 

23 

25 

18 

10 

4 

1948 

420 

262 

351 

378 

360 

340 

322 

316 

277 

268 

260 

351 

15 

20 

26 

23 

1949 

326 

355 

355 

377 

102 


Oxford,  England,  51  ■*46'N,  0I'’l6'W 


Y  ear 

Jan 

Feb 

Mar 

Apr 

Jun 

ZhL 

Aug 

Se£ 

Oct 

Nov 

Dec 

8 

16 

24 

30 

29 

14 

20 

23 

18 

21 

1925 

424 

382 

424 

403 

373 

364 

356 

336 

305 

295 

10 

25 

21 

27 

28 

25 

27 

25 

22 

1926 

333 

347 

358 

375 

367 

318 

309 

273 

277 

25 

29 

30 

24 

24 

25 

10 

17 

1927 

408 

384 

382 

373 

345 

329 

305 

283 

11 

20 

19 

19 

24 

24 

12 

1928 

354 

386 

379 

352 

310 

291 

286 

16 

31 

1950 

303 

320 

31 

28 

31 

30 

31 

30 

31 

31 

29 

31 

30 

31 

1951 

336 

403 

393 

390 

397 

359 

321 

324 

276 

269 

268 

301 

31 

29 

31 

30 

31 

30 

31 

31 

30 

31 

30 

31 

1952 

355 

333 

393 

385 

378 

367 

336 

322 

311 

296 

282 

303 

31 

28 

31 

30 

31 

30 

31 

31 

30 

31 

30 

30 

1963 

303 

345 

355 

379 

355 

364 

340 

307 

283 

279 

268 

313 

30 

28 

31 

30 

31 

30 

31 

31 

30 

31 

30 

26 

1954 

344 

392 

367 

398 

385 

347 

337 

318 

294 

273 

275 

280 

30 

28 

31 

27 

31 

30 

31 

31 

30 

31 

30 

29 

1955 

320 

392 

370 

326 

358 

344 

339 

313 

283 

284 

272 

310 

30 

29 

31 

30 

31 

30 

31 

31 

30 

31 

30 

30 

1956 

350 

369 

373 

389 

362 

358 

341 

340 

294 

286 

280 

286 

29 

25 

31 

30 

31 

30 

1957 

313 

385 

350 

363 

379 

356 

Pasadena, 

Calif.,  34‘09'N,  118 

°10'W 

1939 

306  306  286  286 

1940 

347  354 

354  347  326  279 

Potsdam,  Germany 

1942 

398  373  400  339 

1943 

377  345  379 

279 

103 


Reykjavik, 

Iceland,  64“08'N,  21“ 

54'W 

Year 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul_ 

Sep 

Oct 

Nov 

3 

9 

13 

20 

16 

24 

21 

8 

1 

1952 

341 

423 

409 

398 

370 

345 

330 

325 

288 

0 

16 

13 

30 

23 

28 

21 

21 

10 

1953 

351 

394 

400 

398 

379 

359 

336 

303 

307 

3 

28 

28 

30 

30 

28 

31 

29 

21 

1 

1954 

400 

418 

443 

415 

384 

367 

345 

314 

299 

13 

28 

30 

29 

31 

1955 

398 

392 

394 

377 

354 

1 

Rome,  Italy, 

41”54'N,  12 

“28'E 

1 

1 

1 

30 

31 

30 

31 

31 

30 

31 

29 

31 

1954 

389 

378 

321 

303 

294 

275 

273 

277 

276 

31 

27 

29 

30 

31 

30 

31 

31 

30 

31 

30 

31 

1955 

295 

332 

340 

358 

324 

321 

320 

303 

291 

292 

271 

290 

31 

27 

29 

28 

28 

29 

31 

31 

29 

29 

27 

31 

1956 

326 

382 

348 

360 

348 

366 

325 

311 

292 

295 

305 

316 

30 

28 

31 

26 

31 

30 

1957 

333 

322 

3?0 

362 

366 

332 

Santa  Maria, 

Azores, 

36“55'N,  25 

“lO'W 

27 

3 

1951 

318 

330 

1952 

8 

28 

31 

24 

30 

30 

31 

31 

30 

31 

30 

31 

1953 

294 

311 

345 

356 

343 

341 

306 

299 

279 

271 

275 

283 

29 

28 

31 

30 

30 

30 

29 

26 

30 

22 

30 

30 

1954 

295 

321 

336 

356 

364 

329 

320 

298 

290 

279 

280 

287 

29 

28 

27 

30 

31 

30 

31 

31 

27 

31 

29 

31 

1955 

271 

292 

330 

344 

329 

321 

321 

305 

292 

275 

287 

279 

28 

29 

31 

30 

31 

30 

29 

1956 

313 

322 

328 

351 

360 

332 

317 

104 


Shanghai,  China.  31*11‘N.  121 °2 9 ' E 


Year 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep. 

Oct 

Nov 

Dec 

15 

15 

10 

14 

12 

9 

23 

14 

14 

14 

13 

15 

1932 

288 

330 

356 

339 

328 

326 

307 

311 

317 

307 

301 

298 

6 

4 

a 

9 

'  3 

5 

3  6 

16 

16 

5 

8 

11 

1933 

345 

358 

356 

345 

347 

352 

316 

325 

314 

288 

280 

314 

12 

12 

12 

7 

10 

9 

18 

10 

8 

16 

7 

9 

1934 

341 

337 

364 

339 

337 

326 

317 

306 

301 

294 

268 

296 

12 

6 

12 

4 

12 

5 

9 

3 

10 

12 

9 

10 

1935 

344 

374 

356 

356 

356 

332 

305 

305 

310 

298 

283 

305 

6 

7 

8 

7 

6 

3 

10 

12 

20 

21 

13 

5 

1936 

360 

366 

374 

366 

363 

284 

310 

316 

313 

303 

291 

301 

10 

9 

9 

13 

18 

15 

23 

22 

12 

19 

9 

14 

1937 

318 

343 

351 

362 

360 

335 

318 

302 

302 

296 

302 

322 

11 

10 

15 

11 

6 

16 

22 

15 

16 

24 

15 

1938 

363 

362 

381 

355 

333 

317 

307 

292 

292 

299 

320 

17 

10 

12 

13 

16 

13 

13 

17 

20 

12 

14 

24 

1939 

321 

325 

348 

363 

352 

336 

320 

299 

310 

305 

303 

320 

14 

12 

11 

13 

12 

12 

21 

18 

15 

16 

17 

17 

1940 

348 

364 

384 

385 

360 

363 

320 

324 

310 

303 

302 

318 

21 

13 

17 

18 

15 

17 

16 

16 

10 

20 

16 

15 

1941 

335 

340 

351 

364 

363 

354 

324 

314 

314 

317 

313 

316 

22 

7 

16 

17 

15 

8 

21 

7 

11 

15 

18 

4 

1942 

321 

360 

352 

379 

377 

343 

326 

302 

303 

307 

309 

313 

Srinagar, 

India 

.  34 

•05'N.  74“50'E 

16 

27 

24 

24 

6 

25 

1956 

283 

280 

272 

267 

243 

254 

16 

22 

22 

21 

21 

29 

1957 

275 

264 

268 

275 

282 

273 

105 


Spitzbergen,  Norway,  78 *N _ 15*E 


Year 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

22 

21 

1929 

347 

317 

5 

7 

1950 

310 

306 

5 

2 

6 

13 

8 

8 

12 

8 

1 

4 

1951 

377 

458 

481 

442 

400 

367 

317 

290 

292 

330 

4 

3 

5 

18 

14 

12 

31 

26 

17 

5 

1 

3 

1952 

329 

398 

513 

434 

404 

351 

329 

303 

291 

287 

231 

375 

7 

7 

25 

21 

26 

29 

25 

23 

4 

1 

1953 

404 

396 

407 

428 

388 

328 

302 

277 

231 

400 

4 

11 

31 

25 

30 

25 

25 

31 

27 

26 

3 

1954 

378 

409 

450 

442 

381 

363 

328 

287 

256 

273 

286 

6 

29 

28 

30 

30 

31 

30 

1955 

396 

396 

392 

356 

318 

296 

276 

4 

27 

30 

29 

31 

26 

1956 

408 

434 

386 

347 

301 

284 

Table 

Mountain, 

Calif, 

.  ,  34 

‘'24'N,  117 

'42'W 

3l 

30 

30 

29 

29 

1928 

264 

252 

257 

252 

277 

21 

22 

26 

29 

29 

27 

30  22 

1929 

294 

317 

326 

335 

313 

287 

277  275 

Tateno,  Japan, 

Sb’OS 

'N,  140°08' 

23 

10 

15 

16 

23 

25 

1955 

280 

273 

268 

246 

268 

287 

22 

21 

17 

20 

13 

13 

2 

16 

18 

15 

19 

1956 

322 

348 

337 

348 

320 

311 

314 

287 

273 

261 

258 

6 

29 

30 

30 

29 

1957 

355 

375 

356 

351 

343 

106 


Tromso,  Norway,  69*40'N,  18*57'E 


Year 

Jan 

Feb 

Mar 

Apr 

May 

1935 

9 

7 

17 

25 

28 

1936 

233 

269 

404 

419 

374 

6 

2 

18 

23 

16 

1937 

339 

447 

430 

358 

339 

17 

17 

19 

1938 

394 

426 

382 

21 

24 

21 

1939 

367 

374 

371 

17 

18 

25 

24 

27 

1940 

384 

435 

441 

445 

392 

22 

28 

31 

30 

31 

1941 

352 

443 

469 

443 

412 

18 

28 

31 

30 

31 

1942 

335 

461 

492 

431 

412 

31 

28 

31 

30 

31 

1943 

280 

397 

413 

412 

401 

19 

28 

30 

30 

31 

1944 

180 

379 

403 

400 

382 

31 

28 

31 

30 

31 

1945 

226 

373 

398 

432 

397 

25 

27 

31 

30 

31 

1946 

197 

366 

420 

426 

385 

15 

23 

31 

30 

31 

1947 

167 

384 

409 

404 

389 

10 

26 

31 

30 

31 

1948 

253 

348 

378 

405 

390 

28 

29 

29 

20 

1949 

333 

369 

420 

409 

1950 

13 

14 

10 

1951 

487 

428 

409 

3 

3 

12 

13 

14 

1952 

351 

413 

460 

431 

404 

Jun 

Jul 

Aug 

Se£ 

Oct 

Nov 

Dec 

8 

17 

25 

17 

5 

2 

317 

271 

267 

227 

257 

258 

24 

25 

27 

18 

16 

343 

322 

294 

279 

265 

19 

28 

22 

22 

2 

4 

4 

322 

294 

279 

276 

279 

230 

273 

19 

18 

24 

5 

4 

343 

294 

295 

287 

268 

24 

31 

26 

29 

30 

30 

18 

344 

303 

273 

276 

288 

295 

201 

23 

31 

31 

30 

31 

28 

19 

369 

344 

326 

299 

282 

242 

227 

30 

31 

31 

30 

30 

30 

15 

374 

317 

306 

287 

292 

271 

296 

30 

31 

31 

30 

31 

30 

31 

354 

325 

310 

299 

288 

239 

211 

30 

31 

31 

30 

29 

23 

16 

344 

454 

307 

283 

264 

241 

169 

30 

31 

28 

30 

31 

27 

30 

336 

299 

292 

294 

284 

250 

181 

30 

30 

31 

30 

31 

30 

26 

358 

302 

301 

295 

287 

224 

177 

30 

29 

31 

29 

27 

18 

333 

299 

287 

273 

248 

205 

30 

31 

30 

30 

30 

27 

5 

344 

317 

299 

292 

288 

208 

147 

30 

31 

31 

30 

29 

9 

350 

311 

309 

286 

272 

246 

9 

4 

6 

5 

2 

5 

2 

360 

325 

288 

275 

242 

280 

393 

21 

18 

11 

9 

19 

3 

2 

359 

317 

313 

307 

298 

318 

379 

(Cont.  ) 


107 


Tromso,  Norway  (Cont.  ) 


Year 

Jan  Feb 

Mar 

Apr 

Ma^ 

Jun 

Jul 

Sep 

Oct  Nov  Dec 

28 

31 

30 

29 

30 

29 

30 

1953 

415 

385 

322 

311 

290 

261 

252 

26 

31 

28 

28 

29 

30 

31 

29 

30 

1954 

398 

427 

430 

382 

358 

309 

292 

290 

272 

24 

29 

28 

31 

30 

31 

31 

30 

29 

1955 

456 

437 

416 

393 

364 

313 

294 

268 

250 

26 

27 

29 

29 

25 

1956 

333 

377 

431 

364 

337 

Uccle 

1,  Belgium, 

50° 

48'N, 

04°21 

;e 

7 

25 

13 

21 

1951 

18 

20 

26 

26 

23 

19 

11 

11 

258 

257 

262 

302 

1952 

388 

358 

393 

389 

390 

356 

348 

305 

1953 

14 

17 

19 

19 

18 

18 

14 

9 

16 

6 

1954 

360 

396 

374 

347 

335 

301 

283 

253 

252 

277 

4 

4 

9 

16 

14 

9 

7 

4 

5 

7 

1955 

326 

362 

359 

329 

350 

337 

321 

310 

250 

267 

1956 

7 

1 

3 

1957 

356 

360 

302 

108 


Uppsala,  Sweden.  59‘52'N,  17‘‘37'E 


Year 

Jan 

Feb 

Mar 

Apr 

May 

’  Jtin 

Jul 

Aug 

Sep 

Oct 

Ni^ 

Dec 

23 

16 

10 

6 

1951 

350 

307 

287 

250 

18 

20 

27 

20 

21 

22 

22 

28 

21 

16 

1952 

449 

383 

403 

393 

340 

335 

314 

324 

276 

335 

13 

24 

31 

30 

31 

30 

31 

31 

29 

30 

30 

8 

1953 

309 

373 

364 

385 

398 

358 

359 

330 

287 

279 

269 

282 

24 

28 

30 

30 

30 

29 

26 

21 

26 

30 

23 

1954 

360 

394 

401 

426 

390 

367 

359 

322 

317 

295 

282 

20 

20 

28 

28 

29 

26 

30 

28 

27 

25 

24 

17 

1955 

325 

421 

416 

401 

394 

359 

329 

306 

296 

286 

287 

314 

17 

17 

27 

26 

31 

28 

28 

24 

17 

26 

19 

15 

1956 

360 

351 

379 

415 

384 

358 

339 

341 

279 

288 

290 

292 

19 

21 

27 

24 

25 

17 

1957 

321 

475 

400 

389 

385 

374 

Valentia, 

Ireland,  51 

‘56'N, 

10*15 

I'W 

7 

6 

1926 

273 

303 

19 

16 

20 

12 

14 

13 

8 

9 

1927 

408 

394 

371 

384 

352 

350 

310 

296 

Washington, 

D.  C.  . 

38*51 

'N,  77 

*02' 

W 

4 

6 

no 

12 

17 

2 

1 

18 

15 

12 

17 

1948 

288 

292 

obs 

321 

310 

277 

302 

286 

262 

227 

267 

8 

6 

12 

11 

15 

1 

17 

16 

14 

12 

12 

15 

1949 

276 

284 

318 

330 

325 

374 

295 

290 

276 

249 

271 

305 

9 

7 

5 

1950 

279 

333 

341 

1951 

1952 

1953 

9 

17 

17 

19 

16 

15 

21 

21 

15 

17 

12 

16 

1954 

345 

374 

412 

384 

386 

344 

347 

326 

311 

301 

306 

317 

17 

16 

20 

12 

20 

20 

19 

18 

15 

16 

17 

4 

1955 

3  05 

321 

356 

362 

356 

363 

343 

314 

317 

299 

295 

299 

12 

14 

1 

24 

15 

3 

7 

10 

3 

1956 

337 

277 

281 

371 

329 

322 

326 

320 

268 

